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ABSTRACT 

In  this  paper  the  general  prohlem  of  scattering  of  an  artitrary 
incident  electromagnetic  field  "by  a  conducting  or  non-conducting  obstacle 
is  investigated.  The  partial  differential  equations  and  "boundary  conditions 
satisfied  "by  the  field  are  transformed  into  a  pair  of  inhomogeneous  linear 
integro-differential  equations  for  two  vector  functions.   Por  an  ohstacle 
which  is  a  thin  shell  of  constant  thickness  h,  a  formal  procedure  for  ohtain- 
ing  a  solution  of  these  equations  as  power  series  in  h  is  given.   The  lowest 
order  term  in  this  solution  is  the  incident  field.  An  explicit  expression 
for  the  next  term  is  found  in  the  form  of  a  surface  integral.  This  integral 
is  evaluated  approximately  "by  an  extension  to  surface  integrals  of  the  method 
of  stationary  phase.   The  scattered  field,  given  "by  the  terms  so  obtained,  is 
that  which  is  reflected,  diffracted  or  transmitted  "by  a  thin  shell.  The 
physical  properties  of  this  solution  are  examined  in  detail,  and  satisfactory 
agreement  is  found  with  many  results  previously  o"btained  "by  other  methods. 
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1. 

1,  Introcbiction 

The  problem  of  reflection  and  transmission  of  electromagnetic  waves 
ty  an  obstacle  is  that  of  determining  the  electromagnetic  field  in  a  region  con- 
tedning  a  source  of  radiation  and  an  ohstaxile  with  electrical  properties  differ- 
ent from  those  of  the  surroanding  medium.  The  solution  of  this  problem  can  be  ex- 
pressed as  the  stun  of  two  parts,  the  incident  and  scattered  fields.  The  incident 
field  is  that  field  which  would  be  present  in  the  absence  of  the  obstacle  (i.e. 
obstacle  replaced  by  the  surroanding  medium);  the  scattered  field  is  the  differ- 
ence between  the  total  field  and  the  incident  field,  and  is  due  to  the  presence 
of  the  obstacle.  The  scattered  field  is  elLbo  called  the  reflected,  trsinsmitted, 
or  diffracted  field,  depending  upon  the  region  of  space  in  which  it  is  being  dis- 
cussed. 

Often  the  properties  of  the  source  are  altered  by  the  reflected  field. 
This  effect  (change  in  radiation  impedance  of  the  source)  is  of  importance  when 
the  reflected  field  is  comparable  in  magnitude  to  the  incident  field  at  the  source. 
In  the  present  inrestigation  it  will  not  be  considered,  since  the  reflected  fields 
will  be  much  wesiker  at  the  source  than  the  incident  field.  Fur  the  nn  ore,  it  will  be 
assumed  that  the  region  under  consideration  is  unbounded  and  that  the  medium  is 
homogeneous  and  isotropic,  although  the  methods  axe  applicable  to  more  general  prob- 
lems.  In  addition  it  will  be  assumed  that  the  incident  field  is  known. 

In  order  to  solve  the  problem  it  is  necessary  to  obtain,  in  the  ob- 
stacle as  well  as  in  the  surrounding  medium,  a  solution  of  Maxwell's  equations 
which  satisfies  the  boundary  conditions  on  the  obstacle  boundary.  This  solution 
must  also  contain  the  incident  field.  For  incident  plane  or  spherical  waves  this 
last  condition  can  be  fulfilled  by  requiring  the  solution  to  have  an  appropriate 
singalarity  at  the  source;   in  other  cases  this  condition  may  be  more  troublesome. 

Because  of  the  mathematical  complexity  of  the  problem,  various  methods 
of  approximations  have  been  devised.  Among  these  the  method  of  geometrical  optics 
and  the  Kirchhcff  method  in  the  theory  of  diffraction  are  probably  most  ^ri.dely  used, 
Eecently  R.  K,  Luneberg  (10)  has  shown  (by  the  method  of  characteristics)  that  the 
geometrical  optics  result  is  approached  by  the  solution  of  Maxwell's  ec[aatlons  as 
the  wave  length  tends  to  zero,  for  arbitrary  surfaces  of  discontinuity  in  the  medium, 

» 
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But  long  before  this,  geometrical  optics  was  used  for  tracing  reflected  and  re- 
fracted rays  and  deriving  the  lens  and  mirror  laws.  More  recently  the  relative 
intensity  at  any  point  in  the  field  reflected  from  aJi  arbitrary  curved  surface 
has  been  obtained  exactly  under  the  assumptions  used  in  geometrical  optics  (for 
incident  plane  or  spherical  waves)  (12)(2)  and  the  method  can  be  extended  to  the 
transmitted  field  ejad  other  problems. 

The  geometrical  solution  has  several  limitations.  In  the  first  place, 
it  takes  no  account  of  the  physical  properties  of  the  reflecting  surface  and  there- 
fore can  yield  only  relative  intensities  in  the  reflected  field,  i.e.,  the  ratio 
of  the  reflected  intensity  at  one  point  to  that  at  another.   Secondly,  it  takes  no 
account  of  the  incident  wave  length  and  therefore  does  not  yield  diffraction  effects. 
Thirdly,  it  leads  to  infinite  intensities  at  foci.  In  spite  of  these  objections,  the 
geometrical  solution  yields  useful  Information  about  the  reflected  field,  especially 
for  short  wave  lengths. 

The  second  and  third  objections  mentioned  above  may  be  overcome  by 
using  the  Kirchhnff  method,  which  has  been  widely  used  in  optical  diffrsustion  prob-* 
lems  but  has  been  little  used  in  reflection  problems  until  lately.  It  consists  in 
expressing  the  field  as  an  integral  involving  the  field  and  its  derivatives  on  the 
obstacle  surface,  and  then  assuming  values  for  these  quantities  on  this  surface. 
|(l)  p  yol.Until  now  the  only  justification  for  this  procedure  has  been  that  its 
results  agree  with  experiment.  It  is  believed  that  the  results  of  the  present 
Tjaper  provide  some  theoretical  justification  for  it  in  some  problems  (see  below). 

Rayleigh  (I3)  pp  118-126  proposed  the  use  of  the  Kirchhaff  method  in 
problems  of  au;oustic  reflectioa,  and  several  such  problems  have  been  solved  in  this 
way  by  Willis  (16).  Willis  evaluated  the  integrals  by  making  use  of  the  Fresnel 
zone  concept.  In  the  electromagnetic  case  the  Kirchhaff  method  has  been  applied 
to  reflection  from  sevexal  objects  of  special  shape  at  the  University  of  Pennsylvania 
(17)  .  Finally  the  reflection  of  acoustic  or  electromagnetic  waves  from  an  arbit- 
rary curved  surface  has  been  treated  approximately  by  the  Kirchhoff  method  (12)(7). 
The  solutions  in  the  two  cases  are  identical,  and  yield  for  the  reflected  amplitude, 

In  the  evaluation  of  the  Kirchhoff  integral  in  ref.  (I7)  only  linear  terms  are  re- 
tained in  the  expansion  of  the  phase.  This  is  essentially  the  Frauenhofer  approxi- 
mation of  diffraction  theory.  While  this  facilitates  the  evaluation  of  the  integral 
some  of  the  results  are  physically  ijinreasonable.  The  retention  of  the  terms  of  second 
degree,  in  the  present  p^er  and  in  ref.  (7),  which  is  the  Fresnel  approximation, 
eliminates  these  discrepancies. 
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the  incident  amplitude  multiplied  by  two  faictors.   One  factor  is  the  same  as  the 
reflection  coefficient  obtained  hy  geometrical  optics  (in  directions  of  specular 
reflection)  and  is  called  the  geometrical  factor;   the  other  is  called  the  phase 
factor.  The  phase  factor  approaches  unity  in  directions  of  specular  reflection 
and  zero  otherwise,  as  the  ratio  of  the  ware  length  to  the  obstacle  dimensions 
approaches  zero,  and  generally  has  a  maximum  in  the  specular  directions  for  any 
wave  length.  Thus  an  explicit  demonstration  of  the  law  of  reflection  for  curved 
surfaces  is  obtained.   The  phase  factor  gives  the  resxilt  a  dependence  on  the  ob- 
stacle size  as  well  as  the  wave  length.  At  foci  the  amplitude  remains  finite, 
overcoming  another  defect  of  the  geometrical  result. 

The  Kirchhoff  method  also  hsis  defects.  First,  it  has  no  theoretical 
basis.   Secondly,  it  takes  practically  no  account  of  the  physical  properties  of 
the  obstacle  and  therefore  yields  only  relative  amplitudes  in  the  reflected  field 
and  very  little  about  its  polarization.   It  yields  more  information  about  the  re- 
flected and  diffracted  fields  however,  than  does  the  geometrical  solution. 

In  order  to  overcome  the  defects  of  the  Kirchhoff  solution  it  is 
necessary  to  solve  Maxwell's  equations  inside  as  well  as  outside  the  obstacle. 
Several  special  probleias  have  been  solved  in  which  this  has  been  done.  Probably 
the  first  was  that  of  reflection  of  normally  incident  plane  waves  from  the  plane 
interface  between  two  semi-infinite  media  or  from  an  infinite  plate  of  finite 
thickness  (I5).  Ibe  corresponding  problem  of  acoustics  was  solved  long  ago  for 
any  angle  of  incidence  Iref.  (13 )  pp  86,87  J  ,  but  apparently  this  had  not  been 
done  for  the  plate  of  finite  thickness  in  the  electromagnetic  case  until  recently. 
The  solution  of  this  problem  has  been  given  by  R.  K,  Loneberg,  who  has  also  treated 
the  case  of  several  parallel  jlates  and  of  plates  with  variable  properties  (11 ). 
An  approximate  solution  of  the  acoustical  problem  for  many  parallel  plates  or  for 
plates  of  variable  properties  has  been  given  by  M.  E.  Rose  (lU),  and  an  exact  solu- 
tion for  a  special  case  of  the  plate  with  variable  acoustical  properties  by  Courant 
and  Friedrichs  (5). 

Exact  solutions  have  also  been  obtained  in  the  acoustic  and  electro- 
magnetic cases  for  a  point  source  or  point  dipole  at  the  center  of  a  spherical 
shell  (8)  (9).  Similar  results  can  be  obtained  for  a  line  source  on  the  axis  of 
a  cylindrical  shell. 

The  problem  of  scattering  or  diffraction  of  plane  or  spherical  waves 
incident  on  a  sphere  or  circulau:  cylinder  has  also  been  discussed  by  many  authors 
in  both  the  acoustic  1(13)  PP  236-285  land  electromagnetic  cases  1(15)  p  5631-  The 
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solutions  are  usually  given  in  series  form,  and  these  are  very  slowly  corwergent 
for  wave  lengths  small  comoared  to  the  radius  of  the  sphere  or  cylinder.  If  the 
wave  length  is  large  compared  to  the  radius,  in  the  case  of  a  sphere,  agreement 
with  the  Esyleigh  scattering  law  is  obtained.  The  Eayleigh  scattering  law  gives 
the  field  at  large  distances  scattered  from  an  object  small  compared  to  the  wave 
length  obtained  hy  Hayleigh  for  sound.  The  corresponding  electromagnetic  problem 
has  also  been  treated  1(15)  P.  572]. 

In  the  acoustic  case  the  reflection  of  plane  or  spherical  waves  from 
an  arbitrarily  curved  thin  shell  has  been  solved  approximately  (12),  taking  account 
of  propagation  inside  the  shell,  without  making  assumptions  about  boundary  values, 
but  making  use  of  the  usual  contiraiity  conditions  at  the  boundary.  This  was  done 
by  transforming  the  wave  equation  to  an  integral  equation  which  was  solved  approxi- 
mately. The  solution  for  the  reflected  sound  field  consists  of  the  incident  press- 
ure multiplied  by  three  factors  -  a  geometrical  factor,  a  phase  factor,  and  a  re- 
flectivity factor.  The  first  factor  coincides  with  the  result  of  geometrical 
acoustics  or  optics  in  direction  of  specular  reflection,  the  first  two  fsictors 
coincide  with  the  result  of  the  Kirchhoff  theory,  and  the  third  (reflectivity) 
factor  coincides  with  the  reflection  coefficient  of  a  thin  infinite  plate  except 
for  a  term  involving  the  curvature  of  the  surface. 

In  this  paper,  the  same  method  is  applied  to  the  electromagnetic 
problem,  and  similar  results  sire  obtained.  The  geometrical  and  phase  factors, 
which  are  the  sane  for  all  components,   agree  exactly  with  the  corresponding  factors 
in  the  acoustic  solution.   Since  the  geometrical  and  Kirchhoff  solutions  are  iden- 
tical in  the  acoixstic  axid   electromagnetic  cases  (see  page}©  )  this  agreement  might 
have  been  expected.  Thus  all  properties  of  the  solutions  which  depend  only  upon 
the  first  two  factors  vdll  be  the  same  in  the  acoustic  and  electromagnetic  cases. 
£t  shoxild  first  be  mentioned  tliat  the  solution  provides  a  partial  justification 
for  the  Kirchhoff  theory  as  well  as  for  the  use  of  geometrical  optics  or  acoustics. 
Furthermore,  from  the  geometrical  factor  the  mirror  law  and  conditions  for  the 
formation  of  point  images  or  foci  can  be  obtained,  Tron  the  phase  factor  the  law 
of  reflection  for  curved  surfaces,  the  approach  to  geometrical  optics  in  the  limit 
as  the  wave  length  tends  to  zero  f-r.i   the  change  of  phase  at  a  focus  or  point  image 

can  be  derived.  In  addition,  in  ref.  (12)  the  Kirchhoff  theory  was  applied  to  re- 

the 
fraction  at  the  curved  interface  between  two  media  and  the  results  apply  tq^  electro- 
magnetic as  well  as  the  acoustic  case.  Snell's  law  for  curved  surfaces,  geometrical 
optics  in  the  limit  of  sero  wave  length,  the  lens  law,  and  conditions  for  the  formation 
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of  point    Images   or   foci,    e.s  well   as   the   change   of  phase  at  a  focus,    are  some  of 
the  results   obtainable  from  the  solution. 

The  reflectivity  factors,    on  the  other  hand,   gire  not  the  same  for  all 
components   of  the  electromagnetic  field.      They  are  equal  to  the  reflection  coefficients 
ohtained  for  plane  waves   incident  ohliquely  on  a  thin  infinite  plate  except  for  a 
term  involving  the  curvature  of  the   surface.     Thus,   except  for  this  additionel   term 
conaiderations  of  the  polarization  of  the  reflected  v/ave  would  be   the  same  as   for 
the  Case  of  plane  waves   incident   on  a  thin  infinite  plate. 

The  reflection  of  electromagnetic  waves  from  thick  shells   or  obstacles 
C£in  be  treated  by  the  successive  sroproximation  method  formulated  in  this  paper. 
Since  the  calculations  required  are  quite  elaborate,   however,   an  alternative  method 
of  treating  such  problems,   which  has  been  successfully  employed  in  acoustics, will 
be  described  (see  page  "^G), 

The  method  of  attack  in  this  paper  is  based  upon  a  transformation  of 
the  partial   differential   equations  satisfied  by  the  electromagnetic  field  into  an 
integro-differential   equation.      In  section  2  this  equation  will  be   derived  for  pjn 
arbitrary  obstacle  in  the  field.     We  shall  then  simplify  the  equation  in  section  3 
by  restrictihg  our  considerations  to  an  obstacle  in  the  form  of  a  thin  shell  with 
constant   thickness.     The  remainder  of  the  paper  will  be  devoted  to    the  approximate 
solution  of  the  problem,    in  the  case   of  the  thin  shell,   and  the  discussion  of  this 
solution. 
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2.  Derivation  of  the  Integral  Equation 

We  shall  suppose  that  an  obstacle  of  one  material  is  immersed  in 
an  infinite  medium  of  another  material.  We  shall  designate  the  infinite  medium 
as  well  as  its  material  "bj  (a),  and  likewise,  the  o^bstacle  and  its  material  will 
"both  be  called  (b).  Instead  of  an  ohstacle  in  a  medlian  it  will  be  convenient  to 
think  of  this  situation  as  one  in  which  one  medium  is  embedded  in  another.  Sach 
medium  is  assumed  to  be  homogeneous  and  isotropic  and  may  be  either  a  dielectric 
or  a  conductor,  i.e.,  the  conductivity  can  be  zero  or  have  any  finite  value. 

A  periodic  source  of  electromagnetic  waves,  not  necessarily  a  point 
source,  is  located  in  medium  (a),  suid  the  fields  arising  from  this  source  are 
assumed  to  be  proportional  to  e"   .   (lor  the  treatment  of  other  types  of  time 
dependence  the  method  of  harmonic  analysis  may  be  anployed  since  the  problem  is 
linear.)  Under  these  conditions  it  can  be  shown  |_(15)  P  392 J  that  the/electric 
field  S  satisfies  the  homogeneous  wave  equation 

(  V   *  k  )  E(r)  =  0  outside  the  source  and  the  obstacle,      (l) 

ct 

(  V^  +  k5)  E(r)  =  0  inside  the  obstacle,  (2) 

In  these  equations  k  and  k.  denote  the  values  of  k  in  medium  (a)  and  (b),  respect- 
ively;  k  =  — -    £'  u  ,  where   £'  =  £  +  .  Similar  equations  hold  for  H(r), 

i.e.,  equations  (l)  and  (2)  apply  if  E(r)  is  replaced  by  H(r).  At  the  boundary 

between  (a)  and  (b)  appropriate  continuity  conditions  must  be  satisfied;   these 

conditions  will  be  discussed  later. 

Equations  (l)  and  (2)  can  be  transformed  into  a  single  integro- 

differential  eq^aation  in  the  following  way.  Let  Gr(r,rO  be  a  function  of  r  sm.d  r' 

2    2 
satisfying  the  equation  (V   +  k  )  &(r,r')  =  0,  with  respect  to  the  second  variable 

Si 

r'     thro-jfeliout   soace  except  at  r'   =  r  where  it  has  a  singularity  of  the  form     -; r-. 

ik^jr  -  r'l  1^  -  ^  \ 

Later  we  shall  take  G(r,r')  =  — ; w —  ,  the  Green's  function  for  infinite  space. 

\t  -  tT\ 

Now,  by  Green's  theorem,   if  (J  and  G  are  finite,    single-valued,    differentiable  func- 
tions in  a  three-dimensional  region  V  bounded  by  a  sectionally  smooth  surface  S 
(i.e.,    consisting  of  a  finite  number  of  parts,    each  with  continuously  turning 
tangent  planes),   and  n  represents   the  normal  to  S  into  Y,   then    UU)  p  390  I 


(Q     4t  "  '^  4f  ^      dS       =        /    (Q  1^2  (J  .  5  v2  Q)     dv.  (3) 

This  equation  is  still  vedid  if  Q,  is  a  vector  function,   since  it  holds  for  each 
coinponent  separately.  It   is  understood  that  r'   is  the  variable  of   integration. 
We  now  set  Q,  =  E-E  ,   where  B     represents  the  electric  field  in  the 

ahsence  of  the  obstacle,   i.e.,   the   incident  field,    and  observe  that  E     satisfies 

2  2  ® 

the  equation     (  v      +  k  )  E     =  0  throiaghout   space  except  at  the  source.     We  assume 

8l  O 

that  E  -  E  ,  which  is  the  scattered  field,  has  no  singularities.  Since  E  satisfies 
equations  (l)  and  (2), 

qV^G-GV^<1  =-  k2(E  -  E  )  G  +  k2(E  -  E  )  G  =  0. 
^  a     0      a     0      ' 

if  r'is  in  (a), 

=  -  k2(E  -E)G+k?BG-k2BG 
a     0      D      a  0 

=  (k^  -  kf)  3J. 

if  r'is  in  (b). 

Zrom  equation  (3)  and  the  above  results  it  follows  that,   if  the  region  V  bounded  by 
S  is  in  medium  (a)  and  the  normal  is  directed  into  V, 


^^-v  H-^  ^^^;y)  ^=°.  w 


S" 


Similarly,  if  region  V  bounded  by  S  is  in  medium  (b)  and  the  normal  is  directed 
into  V, 


?  -  ^o>  -H  -  »  ^^^5^)  ^'^^-  ^f)/^iv. 


(5) 


We  now  ^ply  equation  {,k)   to  the  region  V  bounded  by  the  obstacle 
surface  S',  a  sphere  at  infinity  cy   ,  and  a  sphere  C  around  the  singularity 

of  G  at  r'  =  r  if  r  lies  in  medium  (a).   This  choice  of  7  insures  ths.t  neither 

E-E  nor  G  has  a  singularity  in  V.  We  assume,  furthermore,  that  both  E-E 
0  •  0 

and  G  behave  like  outgoing  waves  at  infinity,  and  this  insures  the  vanishing  of 


g. 


the  integral  over  (J      [(l)  p  25  ]  .  The  integral  /  KE-E  )  ~  -  G  ^ — 22)  ^ 

S 

thus  consists  of  two  parts  corresponding  to  <T  and  S'.  The  integral  over  cy 

r  T 


is  ec 


qual  to     -^n(B(r)  -  E  (r)j    (l)  p  24  I  ,  and  hence  we  have  from  (M-), 


E(r)  -  E^(r)  = 


Vf"-»-)^^-^-^^^^^^}- 


(6a) 


if  r  is  in  (a), 


0  = 


if  r  is  in  (h)   .  -^ 


(6b) 


Similarly,   applying  equation  (5)  to  the  obstacle  volume  from  which 
the  point  r  is    excluded  by  a  sphere    <J     if  r  is  in  (b),  we  obtain 


E(r)   -  E^(r) 


s« 


■dS 


2       ,.2 
a_ 


^-^ 


E  G  dV. 


(7a) 


if  r  is  in  (b), 


o^^/|(km)-v-i^^-^  -(-'->:v-') 


dS 


'2       k2 


%-     a 


(7b) 


if  r  is   in  (a) 
The  variable  of  integration  in  each  of  the  integrals  in  equations   (6a)-(7b)   is  r'    , 


9. 


In  equations   (ja)  and  (7b)  the  normal  direction  is   into  (b)  and 
hence  is  opposite  to  that  used  in  (6a)  and  (6h)  which  is  into  (a).     Changing  the 
normal  direction  in  (7a)  and  ijb) ,   and  then  adding,    (6a)  and  (7"b),  and  (6I3)  and 
(7a),  we  have  for  r  in  either   (a)   or   (b),* 

UdS 


c?n                da 

as" 
an 

(s) 

Y 

in  which  E^  and  £     represent  the  values   of  E  on  the    (a)  and  (b)   sides  of  S', 
respectively. 

In  precisely  the  same  way  we  obtain  for  the  magnetic  field 


H(r)  -  E„(r)  =     ^ /lcE»  .  H^  ^  -  «  |S!  - 


(S') 

If  the  point  r  is  on  the  surface  S*  it  is  excluded  by  a  sphere  half 
of  which  lies  in  (a)  and  half  in  (b).   In  this  case  each  of  equations  (6a)-(7'b) 
has  the  left-hand  side  equal  to  —  |E(r)  -  E  (r)  V  ,  When  the  additions  are 

performed,  as  described  above,  we  still  obtain  equation  (8)  (or  (^)^.  Thus  equa^ 
tions  (S)  and  (g')  hold  for  r  anywhere. 

If  the  surface  S'  and  the  volume  T  vanish,  i.e.,  no  obstacle  is 
present,  E(r)  =  E  (r)   (and  H(r)  =  H  (r)).  Thus  the  solution  in  the  absence  of 
the  obstacle  reduces  to  the  incident  field,  as  expected.  The  two  integrals  in 
(S)  (and  (S'))  therefore  give  us  the  reflected,  scattered,  or  disturbed  field. 

In  order  to  evaluate  the  differences  occurring  in  the  integrand  of 
the  surface  integral  in  (8)  (and  (8')^  v;e  may  use  Maxwell's  equations  and  the 
continuity  conditions  that  mast  apply  to  the  fields  at  the  interface  between  the 
two  media.  If  n  and  t  indicate  normal  and  tangential  components,  respectively, 
the  continuity  conditions  are  the  following  \(15)  PP  3^3^,  ^3  J  • 

See  ref.  llsleQuation  (7).   This  equation  is  the  analogue  of  equation  (8)  in  the 

(7) 
present  psf.er ,  except  that/already  includes  houndary  conditions. 


10. 


(9) 


Since 


-»->->       ^     "i;     -^ 

2^=  Sf  -»■  E*^  and  E     =  B°  +  B     ,  we  obtain  from  (9), 

?  _  ?  =  (  i!ji  -  1)  B^  ^,  (11) 

£a  "" 

where  n  is  the  unit  nonoal. 

Similarly,   from  (10), 

i^_?=   (Zll_  1)  H^^.  (12) 

The  value  of  E^-  E  from  (11)  [and  E^-  H  from  (l2)jcan  te  substituted  into  (S) 

and  (8»). 

^B*  SE^ 

The  evaluation  of  the  difference       —  -     ~^-  is  more  difficult  and 

c3n  dn 

requires  the  use  of  Maxwell's  equations.  For- this  purpose  we  set  up  an  orthogonal 
coordinate  system  in  a  neighborhood  of  the  surface  S".  In  this  system  the  coordi- 
aatg  surfaces  are  the  surface  S'  and  its  parallel  surfaces,  and  the  surfaces  gener- 
ated by  the  normals  to  the  lines  of  curvature,      u.    =  const,   and     u^  =  const,    on  S'. 

The  three  coordinates   of  a  point  are  then  the  distance  n  from  S'    along  a  normal 
and  the  values   of     u,    and     vu  corresponding  to  the  intersection  of  the  normal  and 
S'.     The  positive  direction  on  n  i.s  toward  the  medium  (a).      In  the  neighborhood  of 
the  surface  S'   this  coordinate  system  is  orthogonal,   aJid  if  S*    is  given  by  the 
vector      x(    u,,    Up)  having  continuous  partial  derivatives,   it   is  single-valued 

(6)  p  272  I    .   The  so_tiare  of  the  differential   of  arc-length  may  be  written 

?  P  2  2 

ds     =  dn    +  gj^  du  ^     +  g^  d  u^       »  (^3) 

where  g,  and  g_  are  functions  of  n,  u^,  and   Up  \iX3)   P  3^  1  •   I^  i^i  !]_  and  ip 

represent  unit  vectors  along  n, u^ ,  and u  ,  respectively,  the  field  vectors  may 
be  written 

E  =  Ej,  n  +  E3_i^  +  E2i2   ,     H  =  H^  n  +  H^  i^  *  ^2  ^2  *     ^^^^ 


11. 


Since  it  is  assumed  that  there  is  no  free  charge  present,  and 

— i'^t 
the  fields  are  proportional  to  e  ,  Mazvell's  equations  are 


V  •   B  =  0, 


V  •  H  =  0, 


(15) 


SJ  x^=     i^H.         Vx  H  =   (^5LC:  ^  lf£&)Bj  =  _ 


iOL* 


E. 


(16) 


It  is  understood  that  the  differential  operators  in  these  equations  are  expressed 
in  terms  of  the  coordinates  n,  u^,  and  u.  I  (15)  P  ^9  l« 

We  are  now  in  a  position  to  derive  expressions  for  the  differences 
in  the  normal  derivatives  of  the  components  of  S  and  H  as  we  cross  the  surface  S', 
Using  the  first  of  equations  (l4)  with  the  first  of  equations  (15)  we  have  for 
the  electric  vector, 


x/gl    §2 


^  i\^i^^)-  ^  i\/r^)-  3-5^  (^2  TiT)] 


=  0. 


(17) 


The  left-hand  side  of  equation  (I7)  is  the  expression  for  the  divergence  of  E 


in  general  orthogonal  coordinates.      In  order   to  find  the  value   of 


.^n 


aE^ 

d  n 


on  the  surface  we  anploy  equation  (I7)  first  on  the  (h)  side  of  S'  and  then  on 

the  (a)  side  PJid  suttract.  Applying  this  step,  and  carrying  out  the  differentiation 


an 


iE    VSiSp)only,  we  ohtain 


as 


dn 


an 


/i^ 


^n   -^v^  *  af  ^^1^^  ^  'at  ^4  /%) 


'\ 


-   E^^ 


Kf^^h^2'  ai^K^)-  a%(^2Vir; 

It  follows  from  the  first  of  equations  (9)  that  E?  =  E:  and 

a    t 
S^  =  Eg  on  S'.  Hence  equation  (18)  becomes,  at  any  point  (0,u^,u  )  on  S', 


3D.    "    dn 


i^..^=  U^  _  j;t)  ±n^^^ 


n     n     I 


(18) 


12. 


How     -^  f&TEZ    =  2  G    v/gTil         where  &     denotes  the  mean  cuxvature  of  the  sur- 
d  n  ^  T.  2  m  ^    12  m 

face  at   the  point   (0,u^,u  )     I  ref . (3) § IO5J  .     Using  equation  (11)  and  the  fact 


a         h 

that  E.    =  E.    ,   it  follows  that 


SE 


3E 

"an 


-   =    C  ^  -  1)  . 

t-  a 


2G 


(19) 


Qn  an  '    c'  '      n  m 

^  opposite 

(The  mean,  curvature  G     is  positive  if  the  center   of  curvature  is    /    that  side  of 


m 


the  surface  v/hich  agrees  with  the  positive  direction  of  the  normal. ) 

In  eiactly  the  same  way  we  can  obtain  for  the  normal  component   of  the 
magnetic  vector,   using  the  second  of  equations   (I5)  and  the  second  of  equations   (lU) 
together  with  equation  (12), 


SH 


n 


du 


<-if^.. 


an 


)  H" 
/^a 


2&_ 


m 


(20) 


In  order  to  obtain  an  expression  for  the  difference   of  the  normal 
derivatives  of  the  tangential  component  E,    on  the  two  sides  of  the  surface  S', 
we  use  the  first  of  equations    (lb)  with  the  first   of  equations   {ik)  and  equate 
the  u  -components.     For  V  x  E  we  have,   \{13)  P  ^9  j  t 


V2   i 


VxE  =  -^U-E2/i^-     ^\^\ 


+  -i—  1-^  E     -     -^    E„ 
ViTPs    n         an      2 


^h'j= 


y%~ 


^^1^-  ^^^^  '■ 


<3u^ 


] 


n       -2 


Since  Vx  E  -     — ^    H  =  0  on  both  the   (a)  and  (b)   sides  of  the  surface  S',   we 


have  on  equating  Up-components, 


a     f^^ 


a     r,l3 


1  i  w 


^  (^^)  -  d[  K 


2 


/&[ 


-^K^)-  ^^^ 


au^ 


1  COM 

^—^    H^     =  0.      (21) 


13. 


Since,    from  (9)  aad  (10),  eJ  =  B^     and  H^  =  H^  ,   ecpiation  (21)  yields 


(/^.-/^a)    4   *^(a^-    ^)    .       0 


Thus,  while  the  component  H.    is  continuous  across  the  surface  the  normal  dtrivative 

as,  ah 

-r is  not.     Prom  equation  (11)  and  the  fact   that  E.    =  E.    the  last   term  of   (22) 

1                 ^\            ^\ 
hecomes       (l  -     — i —  )       -r •     Hence 


/iT  ^'a  ^-1 


TT   -T5-   -    —  ^rh-/^a^^2   *^(i-7T7)    aH;-    •  (23) 


a 


In  exactly  the  saae  vay,  equating  u, -components,  we  ohtain 


Similarly,  using  the  second  of  equations  (9),  (10),  (ll),  and  (l6), 
respectively,  we  ohtain  for  the  normal  derivatives  of  the  tangential  components  of  H, 

_!5.  .  ^  =  -  iii  (f .  -  ..  )  E*  +  _i_  (i^h.)  JS  (25) 

an   -  3^      c  ^-^h   ^  a^  ^2     ^—    ^^   ^J    au^  •      ^^5) 


1 

an     an      c   -h    a  -1    ^      ^.^       du^ 

Inserting  into  equation  (2)  the  expressions  for  the  differences  that 
appear  in  the  surface  integral,  given  hy  equations  (ll),(l9),(23),  and  (2U),  with 
the  further  assumption  that  the  scales  on  u,  and  Up  are  so  chosen  that  g.  =  gp  =  1 


IK 


on  S',  we  obtsdn  for  the  electric  vector 


]E(r)  =  I^(r)  + 


^^{4!->^^il 


♦  2GG   )m 


^1 


51 


.^r.  -fJ^hK  -  'A>   *  '^  -  -El'  (^iTif  *  "e^' 


1] 


dS 


.^/: 


EG     dV.  (27) 

Similarly,    su'bBtituting  from  equations   (12),    (20),    (25),   and  (26),  we  oljtaln 


E(r)  =  H^(r) 


•'4{<t' 


)     hNII    *     2&G^)n 
n     3  n  m 


•H} 


-i^ 


(C- 


<)   (^1^2 


H&     dV 


1^) 


;i  ./*)(!_  .^t 


A 


•X^. 


au. 


i^4^) 


au 


2  J 


dS 


(28) 


Equations  (27)  and  (28)  are  the  fandamental  integro-differentlal 
eqaations  for  the  fields  in  the  presence  of  an  ohstacle,  and  are  Talld  for  all 
values  of  r.     Fron  these  equations  it  can  "be  seen  that  if  ^^  =  £     and    A.    =   jU 

(in  which  case  k.    =  k   ),   the  integrals  vanish  and  the  fields  reduce  to  the  incident 
fields,   as  expected.     A  formal  procedure,  nemely  a  method  of  iteration,    can  he  used 
for  ohtaining  solutions  of  these  equations.     One  laay  use   the  incident  field  as  a 
first  approximation.     Using  this   first  approximation  in  the  integrals,   a  second 
gpproximation  can  he  ohtained;    etc,     A  modification  of  this  method,  will  be  used 

in  the  following  section  too  obtain  an  approximate  solution  of  the  problem  when  the    • 
obstacle  is   a  thin  shell. 


15. 

3.   Sedaction  of  the  EgogLtions  for  em  0"bBtacle  which  Is  a  Thin  Shell. 

The  integral  equations   (27)  and  (28)  may  te  simplified,   if  the  o^bstacle 
is  a  thin  shell  of  constant  thickness,  "by  transforming  the  volume  integrals  into 
surface  integrals.      In  this  case  we  may  think  of  the  surface  S'   as  consisting  of 
the  two  parallel  faces  S_  and  S__,   neglecting  the  thin  edges  of  the  ohstacle 

which  contrihate  little  to  the  total.     The  surface  integrals  in  eqaations   (27) 
and  (28)  then  extend  over  S^   and  S_j, 

The  integrand  at  a  point  r'     on  S^-  can  he  expressed  "by  a  Taylor's 

series   in  terms  of  its  normal  derivatires  at  the  point  r'     on  S_  iriiich  is  on  the 

normal  passing  throu^  r*    ,  provided  r  is  exterior  to  the  shell.     Thus,   if  A(r' )=A(r,r' ) 

represents  the  integrand  of  the  surfaxse  integral  in  eqaation  (27)  or  (28),    and  n' 
is  normaJ.  into  the  shell  at  r'    ,   on  S^ ,  we  hare 

00 

A(r|i)=r'    ^    ^A(r.)      .  (29) 

m-0  ^n 

in  which  h  is  the  thickness  of  the  shell.  Taking  account  of  the  opposite  directions 
of  the  normals,  on  Sj   and  S--,  since  the  sign  of  /  A(r')  dS  changes  with  the 

direction  of  the  normal,  we  write 

A(r«)dS-H  /   A(r')dS=  /  ^  ^  -^°Mil!l  dS,   (30) 
y,  y,  "iFl^    an'" 

the  negative  tern,  -  /  A(r')  dS,  cancelling  out  the  first  term  in  the  series  (29). 

In  equation  (3O)  the  normal  is  drawn  from  S,  into  the  shell, 

SimilEirly,  if  r  is  exterior  to  the  shell,  we  can  express  the  integrand 
of  the  volume  integral  in  equation  (27)  la  terms  of  its  values  on  S_,  Thas,  with 


16. 


n'  the  unit  inward  normal  to  S.  at  r'  we  have 


/  I(r«)G(r,r«)  dV  =  /  dS  i     ]B(r«  +  >? 'n' )&(r,r«  ♦'^•n')  d  v^' 


^I     ° 


dS 


/"   00  y.i»i  B  r        -1 


'=1    "°-  '-^"  * 


(31) 


c^,^ («->='('.-)  -• 


S  m=0 


The  volume  integral  in  (28)  can  'be  written  in  a  similar  form,  with  B  replaced  "by  H. 
Sulsstituting  from  (30)  and  (3I)  into  (27 )i  with  n  denoting  the  inward- 
drawn  normal  to  S^,,  i.e.,  into  the  ohetacle,  we  have 


.    ,   00   ■u'^'^l    -^ 


^1(^ 

3^  I  s; 


i)(^+  aoG  )  e\ 

'^Q-a.  m'  n 


c- A       Jill       7  an 


+  (k^  -  k^)  G  E^  1  dS. 


2  au. 


(32) 


The  summation  in  the  integrand  of  equation  (32)  steurtw  with  m  =  0,  eren.  though 
that  of  (30)  starts  with  m  =  1,  hecause  the  power  of  h  in  (32)  is  m  +  1.  In 
this  integrand  the  variable  of  integration  is  r'  and  extends  over  S_.  The 
equation  is  valid  only  if  r  is  in  region  (a),  i.e.,  external  to  the  shell.  Sub- 
stituting fro-ii  (30)  into  (28)  and  writing  the  volume  integral  in  the  same  form 
as  (31)  with  E  replaced  hy  H,  we  have 


C& 


^'" = ^o'^'  ^  ^4  S  ^  S 


an  ^  o  '^  an      m   n 


e;  -  £;)  oci/,  -  i^)  *  (1  -  ^"'(Hth^  *  ^2  4? 


au. 


-h 
■^  ) 


+  (k^  -  k^)  G  H^  1  as 


(33) 


17. 

Inspection  of  equations  (32)  auid  (33)  shows  that  if  h  =  0,  E(r)  =  E  (r) 

and  H(r)  =  H  (r),  in  accordance  with  the  fact  that  the  fields  in  the  absence  of 

0 

the  ohstacle  are  S  (r)  and  H  (r).  Since,  on  physicsuL  grounds,  it  is  expected 

that  the  scattered  fields  will  he  small  when  the  thickness  h  is  small,  it  is 
natural  to  seek  fields  valid  for  shells  of  any  thickness  h  in  the  form  of  a  power 
series  in  h.  We  therefore  try  to  find,  in  medium  (a),  solutions  of  the  foim: 

00  00 

E{r)  =  y~  h  B(^)(r),  H(r)  =  ^   h''  H^^^(r)  ,  where  E^^j(r)  and  H(^)(r) 

are  functions  to  he  determined. 

Equations  (32)  and  (33)  can  he  used  for  the  purpose  of  determining  E/  .\ 

and  H/^\«  ^^   order  to  do  this  it  is  necessary  to  transform  the  field  components 

and  their  derivatives  in  the  integrands,  which  refer  to  medium  (h),  into  expressions 

involving  the  corresponding  quantities  in  medium  (a).   This  can  he  accomplished  by 

the  use  of  the  "boundary  conditions  and  Maxwell's  equations,  as  was  done  for  the 

field  components  and  first  derivatives  in  section  2.   If  gifter  this  has  been  done, 

the  expansions  for  E(r)  and  H(r),  valid  in  mediujn  (a),  are  inserted  into  equations 

(32)  and  (33) t  these  equations  are  identities  in  h,  and  hence  like  powers  of  h  can 

be  equated  yielding  equations  for  E/g\  and  H/^w  When  this  is  done  we  find  that, 

m+1 
because  of  the  presence  of  the  factor  h    in  the  integrands,  the  coefficient  of 

h  on  the  right-hand  side  of  (32)  or  (33)  can  involve  only  those  B/  \  and  E,    s  for 

Kl)  vqJ 

which  q  <   v'  ,     Thus  E/^s   and  H,  ,v  can  be  found  as   surface  integrals  involving  these 

E/    \  and  R,    s.     Since  E/    \  =  E     and  H/    s  =  H     are  known,  E/  ,\  and  H/  ,n  can  be  found 

successively  for    ^=  0,1,2,    etc.,   by  integration.     The   expressions   for  Ey^x  and  B./-,\ 

axe  given  below;      for  V  =  2,3,    etc.,    the  expressions  are  similar. 

In  equations  (3^),  (35)  sJid.  (37)   the  field  quantities  with  superscript  b  are 
defined  in  terms  of  the  corresponding  quantities  with   superscript  a  by  equations(9), 
(lO),(l9),(23)-(26),     Therefore  E/    v   and  E,.  may  be  considered  to  be  the  extensions 

of  B/    \  and  H/    \  across   the  boundary  surface  by  means    of  the  boundary  conditions. 

It  is   sufficient  here,  however,   to  consider  them  merely  as  abbreviations   in  order  to 
sinrplify  the  integrands   in  equations    (3U),    (35),    and  (37), 
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«(l/^>^ 


^n  \  V  an  m       n(0) 


^AYa^  ^^44(0r^^l(0))  *  ^^  -  |)  ^^h^'  *  ^2  ^^' 


*     ^<-<^^^}o) 


dS, 


(3^) 


or 


^l/"^=  ^ 


^77  -  1^^ 2       ^m  "ai:^  \(o)  *^—      ^^^^ 

^B.  an  a 


('%-^a'^^^Co)]*h[^Ar«"-''(o)* 


en 


^         f'^^an     az^j^  an  aAj^  ^  B         a'       T.(o) 


+     i. 


m 


-^r^-r.^^^Kio) 


f;-an    a^2 


3A(o)_^ 


ana 


^2 


^^^-^f^'^^ZCo) 


dS 


(35) 


This  expression  for  E/,  x^r)  is  valid  for  r  in  medium  (a).     In  order 

to  evaluate  it  we  must  replace  the  field  quantities  and  their  derivatives  expressed 
"by  means  of  superscript  h     by  quantities  to  which  they  are  equal  in  medium  (a),  "by 
means  of  equations   (9),(lO),(l9),(23)-(26).     Thus,    equation  (35)  "becomes 


19. 


'=(i^'=^/j 


ir,  s(i 


b     an' 


a&> 


^0 


,e"3*^.f5'*(7f-^«i5*=°«.' 


ill*  r»TI* 


S'       ^V 


o) 


*^'-^f^^^'(o)|     *^2 


cTt     /^a''  an  n.(o)  an  c      P)     aj    2(( 


^ 


/      ^V, 


las*/ 


^Hz 


*     (^S  -  ^f)  S  4(0)    »    >     ^ 


(36) 


Similarly,   from  equation  (33)  we  obtain,   for  E.(r), 


H^^-  \. 


^    (A-i)(  ^  + 


20  0   )  H   /    X     n 
an  m       n(o) 


-  ^  U,  -  £^)  aCHi^(„,  -  4  .?(„,)  *  Ci  -^^)  .  (1,  ^  *  i,^  ) 


*     C^-^a)     »=(o) 


dS. 


(37) 
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Carrying  out  the  indicated  differentiation  and  expressing  all  quantities  in  terms  of 
medium  (a), 


H(l/0=     ^/f4^^-^^(4^^^.#     ^^(o) 


dHf 


u  dn  m"'        Bn  mL  ^^  J      n(o) 


/ 


^H)         a'      z^^         n(o) 


+  i. 


itO(  c'        c'U^^  T?^         +  r     '^   2(0) 
-S-^  S  '  ^a^^-d¥  ^2(0)       ^       ^n 


-  ^  ^[A-^>l(o)*^ 


^a 
—  -  1 


as 


3H*. 


'n(o)^   ^  ,  Aa       .  .  _aG      ^^^(o) 


2       Y^b      ■  "^^    ^H 


,U_   ^)a<i-^  *2a    [1.  t^Vtt^  )  Micg  -  k^)  g  Hf,  ,1 


+  i. 


— ^V^a^^a^  ^(0)  *  ^  "air-  *  —  ^  L/'i)-/'ai  ^2(0)*  ^  LeJ  "  ^J       a/(^    ) 


^, 


1    A]  <(o) 


2        ,.  2 


(!<  -  k/)  G     Hj 


2(0) 


dS. 


(33) 


Equations  (36)  and  (38)  represent  a  solution  of  the  problem  of  scatter- 
ing from  a  thin  shell,  correct  to  the  first  powers  of  h.  This  solution  is  valid  at 
all  points  in  medium  (a),  and  is  based  on  the  assumption  that  there  is  no  singularity 
in  the  scattered  field.  Solutions  correct  to  higher  powers  of  h  can  be  obtained  by 
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finding;  expressions  for  l«yp\,  H(p\»  etc.,  "by  a  process  similar  to  that  used  to 

find  the  expressions  for  E/.  v  and  H/.  \  in  equations  (36)  and  (38),  as  outlined  in 

the  paragraphs  immediately  following  equation  (33) •  '^^^  ^^   also  look  for  a  solution 
in  medium  (h)  in  the  form  of  a  series  in  powers  of  h.  The  coefficients  may  he  deter- 
mined hy  a  modification  of  the  method  used  for  finding  those  in  mediujn  (a).  For  the 
investigation  of  reflection,  transmission,  or  scattering,  however,  only  a  solution 
in  medium  (a)  is  necessary,  and  as  we  have  seen,  in  the  case  of  a  thin  shell  this 
can  he  ohtained  without  an   explicit  expression  for  the  solution  in  meditim  (h), 

k.   -Evaluation  of  E/^x(r)  and  H/,  ^(r). 

In  order  to  evaluate  B/,  \(r)  and  H/,  \fr)  the  form  of  the  incident  field 

and  the  expression  for  G  must  he  known.   The  expression  to  he  used  for  G  is  laii  ,,  Ir'-r 
Green's  function  (for  the  wave  equation)  for  infinite  space,  namely  &(r,r')  =  ®_ 


)r.-r| 

which  evidently  satisfies  the  conditions  imposed  in  section  2.  The  expressions  for 

the  incident  field  will  he  written  in  the  form  , 

ik  Ir  -  r  I  ,     ik  Ir  -  r  I 

E(^)(r)  =  d?(r)  e  ^'     °  ,   H^^j(r)  =  h^(r)  e  ^'     °'  .       (39) 

where  r  is  the  center  of  the  radiating  system  that  produces  the  field.  In  particular, 
if  the  source  is  a  dipole  r  is  its  geometric  center;   if  it  is  a  point  dipole  r  is 

its  position;   in  the  case  of  plane  waves  r  is  at  infinity  in  the  direction  from  which 
the  radiation  comes.  It  is  clear  that  any  incident  field  can  he  written  in  the  form 
(39).  We  use  this  form  to  emphasize  explicitly  our  assumption  that  the  main  spatial 
variation  of  the  field  in  the  neighhorhood  of  the  ohstacle  is  expressed  hy  the  ex- 
ponential factor;  that  is,  we  assume  that  c  (r)  and  ^(r)  are  slowly  varying  in  com- 
parison with  this  factor  near  the  ohstacle.  For  plane  waves  (?  and  ^are  constaint  and 
hence  the  asstimption  is  ohviously  correct.  For  other  sources  it  is  the  more  accurate 
the  greater  the  distance  of  the  source.. 

It  will  he  convenient  to  define  the  following  angles: 

•  =  angle  hetween  r'-r  and  n;  thus  cos  Y'=  — '^  ~^Q< 

o  '  an 


P=  angle  hetween  r'-r     and    vn.    ;    thus  cos  (3    =     .^i£ —     ^°  I 

o      1  '        '        sn^ 

^=  angle  hetween  r'-r  and  ML     ;  thus  cos  h     =  3|r'-rQ| 


av<2 


<^-   angle  hetween  r'-r  and  n;  thus  cos  p<  =  9lr'-r  1 

an 
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In  these  definitions  n  represents  the  normal  to  S^   into  the  ohstacle.     With  these 

definitions  and  the  ohservation  that  the  Tariations   in  G    I,    s  and  H/^s  are  mainly 

due  to   the  exponential  factor,   all  the  derivatives  appesxing  in  equations   (36)   and 
(38)  may  "be  approximated  in  the  following  way.     Using  the  chain  rule  we  find 

=   (ik     -   .     ,      .    )  G  coscpc  ^    ik   Gcosc^    ; 
a       I  r '  — r  \  a. 


3G 
an 


=  ik     E,    V   cosr  +     -4^    e 

an  a     Co;  d  n 


ik  lr«-r    ( 

^  °  ^Ik     B     cos   y. 

a     0 


Similarly 


^E 


(0) 


dS 


'^ 


ik    E/    \  cos    (3  ;     -r— ■ 
a     (0;  r  a  u 


(0) 


ik     E/    V   cos    S    ; 
a     (0) 


±±   ^  -k  ^     G  cos^c^   ;     ^ 

2  a  aui  a 


d%(o) 


d  a 


l°n 


-  k       E/    \   cos  ^  cos   &  ; 
a       (0)  '^  • 


2 
a    E 


du. 


122    =  -  k       E/    s   cos  ^  cos   S  . 
a       U; 


(40) 


2     n  ^ 

The  approximations  for   the  derivaties  of  H/    \  axe  similar  to  those  of  E/    s. 

Substituting  from  (39)  and  (UO)   in  equations   (36)jthe  approximate  value 
of  E/^  \(r)  is  found  to  "be 


E(,j(r) 


1 

"57 


r^    6, 


2  2 


(1  -     -S^.)   (  -k       cos   £7<    +  2ik     G     cos  o<) 
eJ  a  am 


+(_^  _  i)(ik  cos^+  2G   )(ik     cosV+  2G„  [l  -  -^1  )  +  (k?  -  kf )  -A- 

+.     [if£(/u^  _  u  )(ik     ^,   COS-*  ik   /f  cos  Y  +  iiifc  '-   £' 1  (f  +  ik^ff  (/^  -  1)   cos  S) 
1  L  c   /  h      ra         aa.  a2  c[d         ajl         an^-^^ 


£h 


_(_^  -  1)  k^^  tf^coso<cos/S  +  (1 ^)(-k^^<?^cos^cos^+2ik^Gjj^  <^n^l-  T^)   ^°^  /^  ^ 
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-  1^  (/<  /x  )(ik    f^   cosc^+  ik   fh.    COS  ^-i±i  ft'  -  (f '1  <?o 

c      /   t      /a  '^al  al  cLt         a-l2 

i-  Ik^  fl-^  [^  -  l]    cos  |6  )  -  (-4  -  1)  \^S^  cosc^cos   S+  (1  -  4)(-k/^  cosYcosS 


+2ik  G    (f    fl  -  4t]  cos  S  )  +(kf  -  k^)    cf, 
a  m     n  L  g  '  J  ^    D         a         i 


ik(lr-r»|  +  lr«-r   |) 


I  r  -  r' 


dS 


(HI) 


The  corresponding  expression  for  H/,  \(r)   can  "be  obtained  from  equation  {hi)  by 

interchanging  <-  with  ^and  £    with  -  U  , 

The   integral   in  equation  (Ul)   is  of  the  type 

f  ifk-  Ir'-rJ    +  k„l  r'-r.l   ) 

Kr^^.r^)  =     /     f(r«)     e     ^  ^  ^  ^  AS 

-4 


(1^2) 


where     f (r' )  = 


/  quantity  in  curly  "brackets  \ 


im  |r«  -  r 


•  V^2  =  ^a'  ^1=^'  ^2=^0- 


and  f(r')  is  slowly  varying  compared  to  the  exponential  factor.  An  approximate  evaluatior 
of  the  integral  in  equation  (U2)  is  given  in  Appendix  I.  In  this  evaluation  the  slowly 
varying  factor  f(r')  is  taken  out  of  the  inte^i^rand  and  given  an  average  value.   Then 
the  integral  of  the  exponential  factor  is  evaluated  hy  a  method  similar  to  that  of 
"stationary  phase",   applied  to  a  double  intet';ral  with  finite  limits  of  integration. 
The  result  of  this  evaluation  is 


K  ^^o"'- v?;'^][^iV^- v^'i^l  • 


(U3) 


The  symbols  in  the  right-hand  member  of  equation  (U3)  have  the  following  meanings 
(see  Appendix  l): 

f  =   average  value  of  f  on  S,. 

/  ^  i  i  H  u^ 

F^(x)  =  /   e    2    du.   In  each  bracket  of  equation  (U3)  the 

signs  agree  and  are  chosen  in  accordance  with  the  rule  given  after  equation  (67  )  of 
Appendix  I. 

^"i  ^  "»  ^  "»  Hi"  ^^   limits  of  integration,  in  suitably  chosen 

units,  corresponding  to  a  rectangular  region  on  the  surface  S  which  accounts  for  the 
main  contribution  to  the  integral. 
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M'  =  Ic.  D,  +  k  D  +  other  terms. 

In  the  ^plication  to  equation  (Ul),  H'  represents  essentially  the  phase  of  the 
reflected  field. 

B.  =  average  value  of  |r'  -  r,  | 

Dp  =  average  value  of  |r'  -  r  | 

In  the  application  these  quantities  represent  avergige  distances  from  the  surface 
to  the  point  of  ohservation  and  the  source,  respectively. 

J  is  given  hy  equation  (65!)  in  Appendix  I,  and  whea  ki  =  kp  =  k  for  the 
application,  it  hecomes 

k   r      2 
J  =  _£:  (  2^)     ^.  _]L  (  1_  [s-sin^r-  sin^c^+  2  8in^«^  sin^  Vsin^  i^-^)] 


-  [cos  cK  +  COS  ^]  [  2G     cos^  "^  ■*■  ^  1 1   (®)  8in^c<]  ) 


2  2  2 

^  (   coBjr-)     ^  J,      (pos^+  cosr)     -  ^  (co8tf<+  COS   V)(2Gj^cos   "*"+&,,  (<^)sin^!r; 


in  which 

0  =  angle  between    u^    and  the  plane  containing  the  normal  and  r^ 

(p=  angle  hetween  u_.    and  the  plane  containing  the  normal  and    r^ 

G     =  mean  curvature  =  •=■(&,   +  G„) 
m  d       J-         c- 

G     =  Gaussian  curvature  =  G,  G_ 
g  12 

G-,Gp  =  principal  curvatures. 

G..    (0)  =  curvature  of  the  normal   section  of   surface  S   deteimined  "by  the  plane 
making  angle  0  with     u, 


2 
(UH) 


2^  .  ,   .,2. 
'2 


=  G.  cos  0  +  G„  sin  9 


The  quantities  in  J  are  determined  at  the  point  of  S  at  which  the  phase 
of  the  integrand  in  ermation  (U2)  is  stationary,  if  such  a  point  exists.  Otherwise, 
they  are  determined  as  described  in  Appendix  I. 

We  apply  the  result  of  the  evaluation  in  eouation  (U3)  to  the  integral  in 
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equation  (Ul)  and  obtain  aB  an  approximate  value  of  E(r): 
E(r)  =  E^Q)(r)  +  Ml^)  r 


=    ^(0)^ 


_¥(n^Ju-4#)( 


2iG 
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l3         a 


-  cos 


c<)  +  (-1  _  i)(i  +  rr^2-— )(i  cos  ^+  ^Fl  -  -is.  I  ) 
£ '  k  cosc3<  k     L  eLJ 
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cos  o( 


cos  Y'L.     1  w 


'Ji^n-A«^^.[^*f3a* 


C08o<J     k  c   cosc^ 
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r£/-,o/Vi£2il^[il^_il) 

Llj     aJ  '         cos«K    n^/Wv         -J' 


*T)  ^a  a  ''  D  k 

a 


+  i. 


a  '  a  L^  J  / 
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fv  V-        2iG  ^  c'     V 

1)  COS  6  cf  -K  (1  -  -11)  ( -  ^2±r+  ^jL^  [i  _  4^iy  cos  s 

'     ^  ^-'n   ^     e'       coso<    k  cos<K  L     £>'  J/^  n 


4 

a 


-  1) 


COSe^ 


ih5) 


The  corresponding  expression  for  H/,  s(r)  can  "be  obtained  from  equation  (U5)  by  inter- 
changing C  with  ^and  £  with  -  u.  . 


Equation  (U5)  and  the  corresponding  equation  for  H(r)  represent  the  electro- 
magnetic fields  in  the  presence  of  the  shell,  and  therefore  give  the  solution  of  the 
problem  within  the  approximations  indicated  in  the  derivation.  In  principle,  higher 
approximations  to  the  fields  could  be  obtained  by  use  of  equations  (32)  and  (33)«  How- 
ever, in  order  to  use  these  equations  to  find  the  second  approximation,  for  example, 
it  is  necessary  to  know  the  fields  obtained  by  the  first  approximation  on  the  obstacle 
surface,  while  the  method  used  here  to  evaluate  the  integral  in  equation  (Ul)  is  valid 
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only  for  otservation  points  at  some  distance  from  the  surface.  Therefore  equation  (U5) 
is  not  valid  on  the  obstacle  surface,  and  cannot  he  used  for  the  determination  of  I/-^ 
and  H/  V.  In  fact,  the  evaluation  of  '£f^\   and  H/^s  on  S  would  require  quite  different 
methods  from  those  we  have  used  and  that  problem  will  not  he  dealt  with  here. 

The  remainder  of  the  present  paper  will  "be  devoted  to  the  interpretation 
and  discussion  of  the  results  obtained. 

5.  Properties  of  the  Solution. 

From  equation  (U5)  and  its  counterpart  for  H(r)  it  can  he  seen  that  the  fields 
as  approximated  therein  are  composed  of  two  parts:  S/  \  and  H,  v,  the  Incident  fields, 

and  hJJ/,  »,  and  hH/,  •>,  the  first  order  scattered,  i.e.,  reflected,  transmitted,  or 

diffracted  fields.  The  scattered  field  components  are  linear  homogeneous  expressions 
in  the  incident  field  coorponenta,  and  hence  if  the  incident  field  is  multiplied  tsy   a 
factor  the  scattered  field  is  multiplied  hy  the  same  factor.  This  was  to  he  expected 
from  the  linear  character  of  the  problem.  We  wish  to  emphasize  the  fact  (which  is 
immediately  evident)  that  each  component  of  the  scattered  field  is  given  as  the  product 
of  three  factors,  each  of  which  has  physical  significance.  The  factors  are  called  the 
geometrical  factor,  the  phase  factor,  and  the  reflectivity  factor,  respectively,  and 
will  now  be  examined  in  detail. 

Geometrical  factor.  The  geometrical  faxitor  G,  the  second  factor  enclosed 
in  large  brackets  in  equation  (U5),  assumes  the  following  form  when  the  expression  for 
J  is  inserted  into  it: 


G  =  <^  1  +  D^ 


f^  fijLl|2_jr    _  tan^^  +  2  tan^c<sin^  Y^sin^  (0  -  ^)] 

\    1   L     cos    << 

2 

-  (coso.  *  COS  r)     [2G^  +  G,,(0)  tan^        *  (^  ^)     *  d/  G^d  * 


cosYv 
cosix' 


^^cosoc^cosr     (2G     cos2Y+G,,(^)  sin2Y)>  "^      _-££l^  (U6) 

\       cos^-  "^  '  J  icos^l 

This  factor  is  the  same  for  all  components  of   the  scattered  fields.      It  is  called  the 
geometrical  factor  for  two  reasons.      In  the  first  place,    it  depends  only  upon  geometrical 
properties  such  as  curvatures  of  the  surface,    distances  of  the  source  aid  observation 
points  from  the  surface,  and  several  angles.      Secondly,   the  amplitude   of  the  field  re- 
flected from  a  surface  has  been  computed  exactly  according  to  geometrical  optics.. 
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(y)   (12),   and  the  ratio  of  the  reflected  aiaplitude  to  that  of  the  field  incident  on 
the  surface  has  heen  found  to  "be  precisely  the  above  quantity  G-  in  directions  of 

specular  reflection  (i.e.  o<=  y',  0  =  (f  )  and  zero  in  other  directions.    (If,    in  the 

2 
geometrical   solution,    one  supposes  only  a  fraction  R     of  the  incident   energy  is  re- 
flected,   then  the  ratio  of  the  amplitude  of  the  reflected  field  to  that  of  the  incident 
field  is  EG.     By  the  use  of   ordinary  geometrical  optics,   however,    one  cannot  determine 
the  value  of  E..  >    This  agreement  leads  to  the  interpretation  of  the  factor  G  in  the 
present  solution  as  that  part  of  the  solution  which  describes  the  geometrical   spreading 
of  the  reflected  wave. 

The   same  factor  G-  eppears  in  the   solution  hy  the  Zirchhoff  method  of  the 
prohlen  of  reflection,   diffraction,   eind  refraction  of  electromagaetic  or  acoustic  waves 
"by  EO-hitrary  curved  surfaces,    (7),(l2).      It  also  appears  in  the  solution  of  the  acoustic 
problem  by  geometrical  optics  and  in  the  more  exact   treatment  of  acoustic  reflection  and 
transmission  by  thin  shells   (12).      It  follows,   therefore,    that  whatever  conclusions  can 
be  drawn  from  the  consideration  of  this  factor  apply  equally  to  all  these  cases. 

In  the  important  case  of  specular  reflection  ( o<  =    ^,  0  =   $2?),  G  reduces  to 
the  form: 

G  =       1  +  2D2   [  5-  -  2Gjjj  cos  ^  -  G       sin  V  tan  •/  J 

or    1  \x  n  -Y  2G- ,,  sinYtan  tT -,  /       ^ 

2     \^2  \  %  D^  J  I        icos^i  •      ^^'^ 

In  this  ca.se  G,   and  therefore  the  reflected  field,   becomes  infinite  v±ien 

5-    =  -       5-    +  2Gj^  cos    )r    +  G  j^    sin  r  tan  V     +     A     ,  (l4«) 


where     A     =  (2Gjj^  cos    Y"  +  G  ^  sin  V^  tan  Y')'^  -  UG  (1^9) 

Points  at  which  G  is   infinite  axe  called  conjugate  points,    and  it  can  be   seen  from 
equation   (HS)   that   for   each  value  of  D,    there  are  two  such  points  along  every  direction 
of  specular  reflection     if  A  is   real  and  not   zero.     When  these  two  coincide   (A  =  0)  a 
point  image  is   fonaed.      The  conjugate  points  and  point   images  are  real   or  virtual  accord- 
ing as  the  sign  of  D_   ,   as   determined  by  equation  (U8)  is  positive  or  negative. 

Point   images   occur   only  if  A  =  0,   and  from  examination  of  the  expression  for 
Awe  see  that   this  can  happen  only  if  (l)  G,    =  Gp  =  0,    i,e. ,  the  surface  point   from 
which  reflection  occurs  is   a  planar  point,    in  which  case  A=  0  for  all    directions; 


It  has   been  shown  by  E.    K.    Luneberg  that   the   reflectivity  factors   of  the  geometrical 
See'ref?  [lO^'p'p'^'yj'jg/'^"  '  ''°^''  '^r.X.rf^,  are  Just   the  well-known  Presnel   formulas. 
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G  1 
or  (2)  G  >  0,  and  if  I  Gj^  1  >  I  Gg  I  ^^^  ®  "  1  •  °08  /  =  (^)^  ,  while  if  |  G^|  <  IG^I 

G,  ^  ■"" 

then  0  =  0  cos  ^  =  (=— )   .  It  is  interesting  to  ohserve  that  the  condition  for  the 

^^2 

occurrence  of  point  images  does  not  depend  upon  the  distance  D,  of  the  source  from  the 

surface  "but  does  depend  upon  its  direction.  In  the  case  of  normal  incidence  cos  <  =  1, 
and  hence  a  point  image  can  occur  only  if  G,  =  Gp,  i.e.,  at  a  spherical  point.  For 

a  point  image  the  image  distance  is   given  "by 

i^  =-  ^  -  2^         .  (50) 

as  can  he  seen  from  equations  (US)  and  (49).  If  the  object  distance  D.  is  infinite, 

the  image  is  called  a  focus,  and  D_  is  then  called  the  focal  length  f. 

Hence 

f=  -^^    .  (51) 

Equations  (50)  and  (5I)  yield  the  mirror  law: 

An  examination  of  the  phase  fsuitor,  to  he  given  later,  will  show  that  usually  a 
change  in  phase  of  180  occurs  at  point  images  and  foci.  A  more  detailed  investiga- 
tion of  the  geometrical  and  phase  factors  can  yield  further  information  ahout  fields 
in  the  neighhorhoods  of  foci. 

The  occurrence  of  such  singularities  as  the  ahove,  vAiere  the  field  hecomes 
infinite,  aiDpears  to  contradict  one  of  the  assumptions  on  which  the  solution  was 
hased,  namely  that  the  reflected  field  is  regular.   The  solution  in  equi^tion  (Ul) 
is  indeed  reg-ilar,  hut  in  the  evaluation  of  the  integral  a  transformation  of  variables 
was  used  requiring  a  division  "by  the  Jacohlan  J,  which  is  essentially  —  ,  At  certain 

points  -  conjiigate  points  and  point  images  -  J  vanishes,  hut  at  such  points  the  trans- 
formation is  not  valid  and  therefore  the  integral  must  he  evaluated  hy  other  means. 
The  present  evaluation  indicates,  however,  that  the  fields  become  large  in  the  neigh- 
horhoods of  such  points.   The  integral  csin  he  evaluated  approximately  at  these  points 
by  other  means,  and  the  results  indicate  that  the  fields  there  are  large  hut  finite. 
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In  the  case   of   specular  reflection  with  normal   incidence   (  <^=     /=  O)  auxd 
D^   =  Dj,,   the   o'baenration  point  is  at  the  source,   and 

G  =     I  (1  +  2G„  D     +  &     D,^)2 
2  m     1  g     1 

1  1 

D     -  2  D       "  2 

=     |(1._1)  (1.     ^)  (53) 

Here  R,    and  E     are  the  principal  radii  of  curvature  at   the  point  on  the  surface  from 

which  reflection  occurs,   and  the  sign  of  G     has  "been  changed  so  that  a  curvature  is 

positive  if  the  center  of  curvature  is  on  the  opposite  side  of  the  surface  from  the 
source.     When  D.    »  IL    add  D,    »  B.     equation  (53)  reduces  to  the  well-known  formula 

G=  i 


2^1  V^ 
Other  special  cases  of  equation  (53)  are: 

Plane  (E^^  =  R^  =  oo  )  :  G  =  p 


1 

1      ^1  ~  2 
Cylinder  (E^  =  oo )  J   &  =  |  (i  +  _£) 

Sphere  (E^^  =  E^)  :     ^  "  |  ^^  *  5"^ 

These  special  cases  show  the  role  of  the  geometrical  factor  in  accounting  for  the 
spreading  of  the  reflected  field, 

Finally,  we  consider  the  case  of  direct  transmission,  («^=  tt  -  ■/,  9  =^). 
Here  -  \ 

«  =  ^p^  (5-) 

This  result  shows  that  G  is  independent  of  the  curvature  of  the  surface.  Moreover, 
the  spreading  acco-onted  for  hy  this  factor  is  the  same  as  for  a  spherical  wave 
diverging  from  the  source. 
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The  Phase  Factor.  The  phase  factor  P,  which  is  the  first  factor  in 
large  "brackets  in  equation  (45),  like  factor  G.  is  the  same  for  all  components 
of  the  first  order  scattered  field.   It  is  giren  "by 


iM' 

P=  ^^ 


This  factor  is  called  the  phase  factor  for  two   reasons.     In  the  first  place,    it 

depends  upon  the  vrave  length  of  the  incident  field  as  well  as  upon  geometrical 

properties,   and  oscillates  as  the  wave  length  is  varied,  with  its  ahsolute  value 

ranging  hetveen  sere  and  approximately  2.      Secondly,    the  problem  of  reflection  and 

refraction  of  electromagnetic  waves  at  arbitrarily  curved  surfaces   and  the  analogous 

acoustic  problem  have  been  solved  approximately  by  the  Kirchhoff  method,    and  these 

solutions  all  contain  the  above  factor  P  as  well  as  the  geometrical   factor  G.     The 

essential  improvement   of  the  Kirchhoff  method  over  that   of  geometrical  optics   lies 

in  the  fact   that   the  former  takes  account, of  the  phase  while  the  latter   does  not. 

ika  (Di  +  EgO 
In  addition,  P  contains  the  factor     e  which  describes   the  variation  in 

the  phase  of  the  scattered  wave  due  to  the  total   distance  traveled.     For  these 
reasons  the  factor  P  introduced  by  the  Kirchhoff  theory  may  be  called  the  phase 
factor.     This  factor  also  occurs  in  the  more  general  treatment  of  acoustic  reflection 
and  transmission  Isee  ref.    (12) J  • 

As  is  shown  at  the  end  of  Appendix  I,   for  a  fixed  direction  of  the  inci- 
dent field  the    ihase  factor  has  its  maximam  absolute  value   in  the   direction  of 
specular  reflection  for  an  obstacle  that  is  symmetrical  about  the  point   of  re- 
flection.     This  is  in  accord  with  the  experimental  result  that  the   reflected  field 
is  stronger  in  the  specular  direction  than  in  non-specular  directions,   and  indicates 
that  the  phase  factor  accounts  for  this  difference.      In  ref.    (12)   it  is  shown  that 
a  slightly  modified  phase  factor  occurs   in  the   solution  of  the  refraction  problem 
stnd  has  its  maximum  in  the  direction  determined  by  Snell's  law.     As  the  wave  length 
tends  to   zero   (k     ->  oo  ) ,    the  phase  factor  approaches   zero  in  all    directions  except 

SL 

these,    i.e.,    the  direction  of  specular  reflection  in  the  csise  of  reflection  and  the 
direction  of  Snell's  law  in  the  case   of  refraction,    for  which  its  absolute  value 
E^proaches  unity.      This    is     shown  in  Appendix  I.      Th'ls,   in  the  limit,   oar  result 
yields  the  lavf  of  reflection  which  is  usually  taken  as  an  assumption  in  geometrical 
optics   (it  is  sometimes   deduced  from  the  assu;nption  of  the  principle  of  least   time). 
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The  law  of  reflection  has  often  bean  decjaced  from  electromagnetic   theory  for  the 

reflection  of  plane  waves   from  an  infinite  plane;      the  present   derivation  also 

applies   to   ciirved  surfaces  and  is  not   restricted  to  plane  waves.      Similarly, 

the  modified  phase  factor,  which  is   obtained  when  the  Kirchhoff  theory  is  applied 

to  refraction,  yields  Snell's  law  for  curved  surfaces. 

In  Appendix  I    it  is   shown  that  the  symbol  F  (x)   is  written  F. (x)  or  F  (x) 

—  p  ~ 

according  as  A  is  positive  or  negative  \Aere  A  is  the  coefficient    of  x     in  the  phase 

of  the   integrand  of  actuation  (62).      Similarly  the   sign  in  F   (y)   is   the  same  as   that 

2  ~ 

of  C^    the  coefficient  of  y     in  the  phase  of  this  integrand.     By  examining  the   ex- 
pressions for  A  and  C  given  in  Appendix  I   it  can  be   seen  that  on  the  specularly 
reflected  ray  through  a  point   image,    both  A  ahd  C  have  the  same  sign  on.  one  side 
and  both  change  sign  on  crossing  the  point  imsige. 
Thus   the   expression  for  the  scattered  field  in  equation  (U5)   contains   the  Fresnel 

integrals  F^  on  one  side  and  F__  on  the  other   side  of  the  point  image.      If  the  argu- 

"                                     1       -''S' 
ments   of  these  functions   are  large,    then     F    ^ e        .     Thus  the  nroduc.t  of 

,  +  i  n        r2      _i  TT 
12  12 

Fresnel  integrals  changes  from  ^  e        to       —  e      on  crossing  a  point 

image  silong  the  reflected  ray.  Thus  there  is  a  change  of  phase  of  180°  in  the 
scattered  field  on  passing  throu^  any  point  image,  in  particular  a  focus,  formed 
by  reflection.   It  is  shown  in  ref .  (12)  that  the  same  considerations  ^ply  to 
point  images  fanned  by  refraction. 

The  arguments  of  the  Fresnel  integrals  are  given  by  the  following  ex- 
pressions; 

^^  and  yi^   are  given  by  similar  e:q)ressions  in  which  ^'  and  ri  '  are  replaced, 
by  ^'  and  >^  '  ,  respectively.   In  these  expressions  ^'  ,  ^,',  >7'  ,  Vi '  are 

limits  of  integration  over  the  rectangular  region  eouivalent  to  the  surface  and  re- 
ferred to  the  point  of  specular  reflection  as  origin.   The  denominators  in  these  ex- 
pressions can  be  given  a  simple  interpretation  in  terms  of  the  Fresnel  zone  concept, 
as  follows. 
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Radiation  from  all  points   in  a  certain  neighborhood  of  a  given  point 
on  the  surface  arrives  at  a  point  r    with  phase  differences  lesB  than.    fl.      There- 
fore the  radip.tion  from  any  point  of  this  neighhorhood  reinforces   that  from  other 
points   of   the  neighhorhood,   rather  than  cancelling  it.     Such  a  neighborhood  is 
called  the  first  Fresnel  zone  around  the  given  point.     If  the  radiation  comes 
from  the  point  r       and  r'    is  a  siirface  point,   then  the  phase  of  the  radiation 

reaching  r  froii  r'   is  Ic.    I  r'   -  r|     *  K     I  r»  -  r    |   ,   iidiile  the  phase  of  that  from 

the  given  point  is  k-D,    +  ^^p*     ^o'^  points  on  the  "boundary  of  the  first  Fresnel 

zone  the  difference  of  these  vaJLues   is     *  TT  .     Using   the  coordinate  system  given 
in  Appendix  I  with  origin  at  the  given  point  on  the  surface  and  the  expressions 

for  lc,|r*-rl     ^k,     jr'-rl     given  in  the  exponent  of  the  integrand  in  equation 

(62),  we  have  as  the  equation  of  the  hoxindary  of  the  first  Fresnel  zone, 

Ax^  +  Bxy  +  Oy^  +  Dx  +  I&r  =     ^  tt     .  (57) 

Rotating  and  translating  the  coordinate  system 

A'x'^     f     C'y'^     =     i    TT     , 
or 

I    A '  I  I  C  M 

In  equation  (58)  the  sign  is  .chosen  the  same  as  that  of  A'  C  .  Thus  we  see  from 
equation  (58)  that  the  boundary  of  the  first  Fresnel  zone  is  an  ellipse  if  A'C  >  0, 
a  hyperbola  if  A'C  <  0,  and  two  parallel  lines  if  either  A'  or  C*  is  zero.  From 
equation  (58)  it  is  seen  that  ^/^Yi  ^^^  ^/ffTT  ^^^   *^®  half-dimensions  of  the  first 

Fresnel  zone  around  the  given  point.  Thus  the  arguments  of  the  Fresnel  integrals 
are,  except  for  the  factor  ^fc^,   the  ratios  of  the  half -dimensions  of  the  surface 
to  corresponding  half-dimensions  of  the  first  Fresnel  zone. 

When  the  dimensions  of  the  surface  are  large  compared  to  those  of  the 

first  Fresnel  zone  the  arguments  of  the  Fresnel  integrals  are  leurge,  and  the  integrals 

1    ±  l^^ 
may  be  replaced  by  their  asymptotic  values,  e    ^   .   In  this  case  the  phase  factor 

r2 

becomes  i  eiM'  multiplied  by±l  or  -  I  ,  the  choise  of  multiplier  depending  upon  the 
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signs  in  the  Fresnel  integrals.  From  the  form  of  M«  it  can  then  "be  seen  that  if 
mediTim  (a)  has  zero  conductivity,  and  k  is  consequently  real,  M'  is  real  and 

EL 

therefor^  the  phase  factor  has  its  absolute  value  equal  to  one.      If  the   surface 
is  much  larger  than  the  first  Fresnel   zone  it  will  contain  many  Fresnel  zones. 
This  is  usually  the  case  unless  the  surface  dljnensions  are  small  compared  to  the 
wave  length,   or  one  or  "both  of  the  radii  of  curvature  are  infinite,   i.e^   the  sur- 
face is  cylindrical  or  plane. 

When  the   surface  is  small   compared  to  the  first  Fresnel   zone  the  argu- 
ments  of  the  Fresnel  integrals  are  small,   and  these  integrals  may  he  approximated 
by  the  first  terms  in  their  Taylor  expansions,  which  yields  ?. (x)  %  x.     Then 


From  equation  (59)  we  see  that,   apart  from  the  exponential  factor,   P  is  equal  to  the 

quotient  of  area  of  the  obstacle  surface  divided  by  the  approximate  area  of  the 

first  Fresnel  zone. 

Since  both  A'   and  C  have  the  faisrtor  k   ,  which  is  equal  to     -^  ,  where 

^  '^a 

X     is   the  wave  length  in  medium  (a),    if  the  conductivity  of  medium  (a)  is  zero,     P 

is  inversely  proportional  to  A     .     The  reflectivity  factor  also  has  the  factor  k  h. 

a  a 

and  hence  one  sees  that  the  scattered  field  is  proportional  to  the  obstacle 
volume  (=  area  xh)  and  inversely  proportional  to  X  ,  This  is  the  analogue  of  the 
Rayleigh  Scattering  Law  for  electromagnetic  waves.   In  particular,  it  is  valid  if 
X  »  obstacle  dimensions,  since  then  the  first  Fresnel  zone  is  larger  than  the 

obstacle. 
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Eeflectivlty  Factor.     The  scattered  field  is  eacpressed  in  eqaation  (U^) 
in  terms  of  conrponents  in  the  n,   i-^,   ig  directions,  where  n,   i^^,   ig  are  nmtually 

orthogonal  unit  vectors   at   the  point  of  stationary  phase  on  the  surface  S.     Each 
of  these  components  is  the  product  of  the  geometrical  and  phase  factors  and  a  third 
factor  which  is  called  the  reflectivity  factor.      Since  this  factor  is  not  the  same 
for  all  components,   it  is   appropriate  to  speak  of  three  reflectivity  factors,    one 
for  each  component.       The  factors  are  the  coefficients  of  n,    i^,    and  i^  in  equation 

(1+5). 

The  reflectivity  factors  are  linear  and  homogeneous  in  the  components 
of  the  incident  field.   They  depend  upon  the  physical  properties  of  "both  media, 
the  thickaess  and  curvature  of  the  shell,  and  the  angles  of  incidence  and  reflection. 
Thus  it  can  he  seen  that  these  factors  take  account  of  the  propagation  within  the 
shell.  From  an  examination  of  the  form  of  these  factors  it  is  seen  that  they  sa"e 
proportional  to  the  amplitude  of  the  incident  field,  and  increase  in  ahsolute 
value  as  U '  -  €i\   and/or  lu  ~ /*y.\    increase.  Moreover,  they  are  proportional 

to  k  h  which,  in  the  case  when  medium  (a)  has  zero  conductivity,  is  -r —  .  Thus 
a  Aa 

the  scattered  field  increases  with  the  ratio  of  the  thickness  to  the  wave  length. 
This  inverse  dependence  upon  X   ,   in  conjunction  with  the  phase  factor,  leads  to  the 

electromagnetic  analogue  of  the  Hayleigh  Scattering  Law,  as  has  "been  pointed  out 
ahove.  When  reflection  and  diffraction  prohlems  are  treated  hy  the  Kirchhoff  method 
no  reflectivity  factor  is  ohtained,  and  hence  the  Rayleigh  Scattering  Law  cannot  he 
deduced. 

If  the  reflectivity  factors  are  divided  respectively  hy  the  corresponding 
field  components,  they  agree  with  the  exact  coefficients  obtained  "by  E.  K,  Luneherg 
(13)  for  the  reflection  of  plane  waves  from  a  thin,  infinite,  plane  plate,  except 
for  additional  terms  proportionfJ.  to  the  mean  c\xrvature  of  the  surface.  The  polariza- 
tion effects  are  therefore  prgictically  the  same  as  for  the  plane  case  ahove.   If  the 
result  in  equation  (U5)  is  specialized  to  the  case  of  plane  waves  incident  on  a  plane 
shell,  the  geometrical  factor  and  the  phase  factor  hoth  hecome  unity,  and  the  re- 
flected and  transmitted  fields  agree  exactly  with  those  obtained  "by  Luneherg,  This 
is  shown  in  Appendix  II  . 

The  presence  of  the  quantity  cos  ^  in  the  denominator  of  some  terms  of 
the  reflectivity  factors  appears  to  indicate  that  the  scattered  field  "becomes  infinite 
as  o<  ->  •p  •  In  this  case,  however,  the  method  used  to  evaluate  the  integral  is  not 
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applicable.      It  should  "be  noted  that  Luneberg's  axact  resiilt,  wMch  remains 
finite  as  o^->  -s  ,   also  contains  cos  '^^  in  a  denominator  when  it  is   expanded 

to  apply  to  thin  shells.      This  shows   that  the  expansion  is  not  valid  as<x  ->  ^  , 

just  as  the  method  we  have  used  for  evaluating  the  integral  fails    in  the  same 

case. 

6.  Conclusion. 

In  this  paper  the  general  prohlem  of  scattering  of  an  arbitraxy  inci- 
dent electromagnetic  field  hy  an  obstacle  has  heen  transformed  to  the  solution  of 
a  pair  of  inhomogeneous  linear  integro-differential  equations.  For  an  obstacle 
which  is  a  thin  shell  of  constant  thickness  h  we  have  given  a  formal  procedare 
for  obtaining  solutions  of  these  equations  as  power  series  in  h.   The  lowest  order 
term  in  this  solution  is  the  incident  field,  sind  an  explicit  expression  for  the 
next  term  has  been  found  in  the  for:a  of  a  surface  integral.  This  integral  has 
been  eveiluated  approximately  by  an  extension  to  surface  integrals  of  the  method 
of  stationary  phase.   The  resulting  expressions  for  the  first  order  scattered 
"ield  give  the  field  reflected,  diffracted,  or  transmitted  by  a  thin  shell. 

This  solution  has  been  compared  with  many  res-alts  previously  obtained 
by  other  methods,  and  in  all  cases  the  agreement  has  been  satisfactory.  Since 
the  previous  results  have  been  extensively  checked  by  experiment,  the  agreement 
furnishes  some  justification  for  our  approximations.  On  the  other  hand,  our  result 
has  been  obtained  directly  from  the  fundamental  equations  of  electromagnetic  theory. 
The  above  comparisons  therefore  ahow  that  the  results  of  the  other  methods  are  con- 
seqaencss  of  these  equations,  and  also  show  under  what  conditions  they  are  valid. 
Most  of  these  comparisons,  such  as  those  with  results  of  geometrical  optics,  the 
Kirchhoff  nethod,  and  the  Eayleigh  Scattering  Law,  were  discussed  in  section  5, 
Some  of  these  results  have  not  heretofore  been  deduced  from  electromagnetic  theory. 
The  solution  leads  to  a  rather  complete  qualitative  as  well  as  quantita- 
tive understanding  of  the  phenomena  of  reflection  and  transmission  by  thin  shells. 
It  also  gives  some  insight  into  the  relationships  among  the  various  approximation 
methods  which  liave  been  used  to  treat  these  problems..  Thus  one  finds  that  in 
geometrical  optics  the  field  reflected  to  any  point  comes  entirely  from  one  point, 
the  stationary  phase  point,  on  the  surface  of  the  reflecting  object.  Since  making 
the  phase  stationary  is  equivalent  to  making  the  optical  path  stationary,  one  feees 
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here  the  reasoa  for  the  appearance  of  Fennat's  "principle  of  least  time"  in. 
geometrical  optics.  In  the  Kirchhoff  theory  the  field  specularly  reflected  to 
any  point  comes  mainly  from  a  neighborhood  of  the  stationary  point  on  the  surface 
of  the  reflecting  o'oject.  Tnis  accounts  for  the  utility  of  the  Tresnel  zone  con- 
cept in  such  proTslems.   Finally,  in  the  present  paper,  the  main  contriTxition  to 
the  field  specularly  reflected  or  transmitted  to  any  point  comes  from  a  volume 
neighhorhood  of  the  stationary  poiat^ inside  the  shell,  as  well  as  from  a  neighbor- 
hood on  the  surface  of  the  shell.  This  is  the  reason  for  the  appearance,  in  the 
present  solution,  of  the  reflectivity  factor  vdaich  takes  account  of  the  electro- 
magnetic properties  of  the  shell  material. 

For  reflection  from  thick  shells  additional  terms  in  the  series  solu- 
tions would  have  to  Tse  evaluated.  However,  the  necessary  calculations  "become  ex- 
tremely complicated  as  was  explained  previously  (see,  also,  comment  in  the  next 
paragraph).  An  alternative  roethod  of  dealing  with  this  problem  practically,  which 
was  found  useful  in  acoustics,  is  the  following.  From  the  physical  significance 
of  the  three  factors  in  the  solution  of  scattering  "by  a  thin  shell  it  may  he  ex- 
pected that  the  field  scattered  by  any  obstacle  will  also  contain  three  such 
factors.  This  is  farther  indicated  by  the  fact  that  the  phase  and  geometrical 
factors  can  also  be  obtained  by  the  Kirchhoff  method  and  the  geometrical  factor 
by  geometrical  optics,  and  both  of  these  solutions  apply  to  any  obstacle  -  not  only 
to  a  thin  shell.  Therefore  one  might  expect  that  these  two  factors  will  have  the 
forms  already  fo\ind  for  them,  depending  only  upon  the  "front"  surface  of  the  obstacle. 
The  reflectivity  factor  however,  depending  upon  the  thickness  as  well  as  the  shape 
of  the  rear  surface,  may  have  quite  a  different  form  from  that  which  it  was  found  to 
nave  for  thin  shells.   Since  in  this  case  the  reflectivity  factor  is  practically  the 
same  as  the  reflection  coefficient  given  by  R.  K.  Luneberg  for  the  reflection  of 
plane  x^aves  from  a  thin  infinite  plate,  it  is  not  unreasonable  to  assume  that  in  the 
case  of  a  thick  shell  it  will  be  given  by  the  corresponding  expression  for  reflection 
of  plane  waves  from  a  thick  plate.  Calculations  based  on  these  assumptions  seem  to 
agree  with  experiment  in  the  acoustic  case  (12), 

In  order  to  provide  a  farther  test  of  the  result  of  this  paper,  an 
exact  solution  for  the  reflection  and  transmission  of  spherical  waves  by  a  spherical 
shell  with  a  dipole  source  at  its  center  has  been  obtained  explicitly  (9).   (it  is 
remarkable  that  the  solution  of  a  problem  apparently  so  simple  leads  to  so  compli- 
cated a  result).  This  solution  was  soecialized  to  the  case  of  a  thin  shell  and 
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compared  to  the  result  of  the  present  paper  vdien  applied  to  this  case.  The  agree- 
ment 'betv/een  the  two  solutions  was  satisfq.ctory  for  all  components  except  E  ,  the 

radial  component  of  electric  field,  and  except  near  the  sorarce.   It  is  at  first 
wight  surprising  that  any  agreeTaent  at  all  was  obtained,  since  in  this  case  the 
reflected  field  has  a  singularity  at  the  origin  and  therefore  the  result  of  this 
paper  does  not  apply.  The  agreement  may  "be  understood  by  assuming  that  the  result 
of  this  paper  was  applied  to  the  problem  with  a  slightly  eccentric  source,  in  which 
the  reflected  field  does  not  have  a  singularity;   the  solution,  except  near  the 
source,  should  be  the  same  as  that  obtained  with  the  source  at  the  center.  The 
error  in  the  small  longitudinal  component  S  is  due  to  the  approximate  evaluation 

of  the  integral  in  equation  (Ul),  to  which  E  is  very  sensitive  in  this  case.  On 

the  whole,  the  agreement  is  quite  encouraging,  indicating  that  the  results  of  this 

paper  may  even  be  valid  in  cases  when  there  are  singolarities  in  the  scattered  field, 

except  near  these  singularities.  The  sensitive  dependence  of  E  is  unimportant, 

longitudinal         ^ 
since  in  most  electromagnetic  fields  the/      components  are  negligible. 
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Appendix  I. 

In  this  appentLix  we  shall  ottain  an  approximate  evaluation  of  the  integrgj. 

,  i(]c,    |r«-r-l      +  k_    |r«-r-|    ) 

Kr^.r^)  =     jg  f(r')  e      ^  ^  ^  ^  dS»      .  (U2) 

where  r'    is  a  point  on  the  surface  S  and  r, .r^  are  any  two  points  many  wavelengths 

away  from  the  surface.      This   is  the   type  of  integral   that   it  is   necessary  to   evalu- 
ate in  order  to  ohtain  the  first  order  scattered  field  in  equation   (Ul).      For   the 
purposes   of  the  present  paper  it  would  "be   satisfactory  to  have  k     =  k   ,   and  this 

specialization  will  he  made  when  the  result  of  this  appendix  is  gpplied  to  the 
integral   in  equqtion  (Ul).     Hov.-ever,    since  an  integral  of  the  form  given  in  equa- 
tion (1+2)  with  k,    ^  k     occurs   in  the  treatment   of  refraction  "by  the  Kirchhoff  method, 

it  is   convenient  to  retain  k,    and  kp,   and  therehy  deal  with  this  case  also.     It   is 

assumed  that   f(r' )   is  slowly  varying  compared  to  the  exponential  factor.     This   is   in 
accordance  with  ohservations  ahout   the  expression  for  f(r')   in  equation  (U-l). 

We  construct  a  rectangular  coordinate  system  with  its  origin  at  any  point 
on  the  surface  S,   having  the  z-axis  nonnal  to  the  surface  and  the  i.y-axes  in 
directions   of  principal  curvature  of  the  surface  at  the   origin.     Vith  this  choice 
of  coordinates,   the   equation  of  the  surface  is 

z  =  ax^  -^  by^  + (60) 

with  no  term  in  xy.     If  D,   and  D     are  the   respective  distances  from  r,    and  r^  to 

the  origin,    one  may  expand     r'  -r,      for  r'   near  the  origin,    as  follows: 

1 


|r.  -  rj   =   [(x-x^)2My-^)2  .  (..,^)2]  2  ^  p^  _     ^  -  -^ 


■d;;!  2 --1-2^57^   r      -^1 


^     'i-az,    -i(^)     U^-^ 


2  -  ^^-   2^5;^^ 


2     Vl 
y =-  xy 


+  terms   of  higher  order.  ((i-\\ 

In  this  expression  x- ,   y.,    z^  are  the  coordinates  of  r^ .     A  similar  expression 
may  he  obtained  for    |r*   -  r_|  hy  replacing  the  subscript  1  by  the   subscript  2  in 
equation  (6I).     In  the  expression  so   found  x  ,  y       z„  are  the  coordinates  of  r_  . 
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If  terms  up  to  the  second  degree  in  the  expansions  for|r'  -  r.l 

andjr'  -  r„l are  inserted  into  equation  (U2),  one  obtains  the  following  approxi- 
mate expressions  for  l(r, ,r  ): 

Kr^.r^)*^  f(r.)e^  t^'  *  ^^'  -h  Cy^  .  Dx  *  ^  .  m]  ^ 

In  this  equation  the  following  abbreyiations  have  been  used: 

2  -,  ,        .  _  2 


(62) 


\ 


'=\ 


[I  -  ^^1  -  1^72^  1    *     d|    [I  -  ^^2  -  1^72^] 


-i 


1  2 

2  n  ,.       r  „  2 


B  =  - 


^Vl         ^2^2^2 


.,3  .^3 


^"-     ^--l^-  =^-^-^-       "='^■'1*^2 

From  equation  (62)  it  can  be  seen  that  the  phase  of  the   integrand 
i.e.    the   exponent,  will  be  stationary  at  the  origin  if,   and  only  If,  both  D  and 
E  are  zero.      If  there  is    a  point    on  S  for  which  both  D  and  E  are   zero  when  that 
point   is  chosen  as  origin  of  coordinates,    then  the  nei^borhood  of  that  point  will 
give  the  main  contribution  to   the  integral.     This  is   the  principle  of  stationary 
phase,   and  is  the  more  valid  the  larger  are  k^  and  kp,  which  in  the  present  case 
meeins   the   shorter  are  the  wavelengths     X,    and  \„.     The  condition  D  =  E  =  0  is 

equivalent  to  the  law  of  reflection  when  k^  =  k_,   and  to  Snell's  law  when  k,    f  kp. 

This  meajis,    in  the  case   of  reflection  of  radiation  from  a  source  at  r     to  a  point 

r       that   the  main  contribution  comes   from  the  neighborhood  of  that  point  r«   on  the 

surface,    if  any  such  point  exists,    for  which  the  normal  to  S  lies   in  the  plane  of 
r,  ,   rp,   r'   and  bisects  the  angle  between  r'  -  r^^  and  r»  -  r   .     In  the  case  of  trans- 
mission by  a  shell,   we  have  k,   =  kp,    and  the  condition  D  =  E  =  0  means  that  the  main 
contribution  to  the  wave  trginsmitted  to   r,    from  the  source  at  r     comes  from  the 
neighborhood  of  the  point  in  which  the  straight  line  from  r.    to  rp  cuts  the  surface, 
if  there   is   sxich  a  point. 
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Whenever  such  a  stationary  point  can  he  found  on  the  surface,  it  is 
to  "be  chosen  as  the  origin  of  coordinates.  In  case  there  is  no  such  point,  the 
principle  of  stationary  phase  indicates  that  the  integral  will  be  smaller  in 
absolute  value  than  it  would  he  if  there  were  one;   this  is  home  out  "by  the 
ohservetion  that  non-specular  reflection  is  much  weaker  than  specular  reflection, 
and  will  he  explicitly  shown  hy  our  result.  In  this  case  a  point  for  which  the 
phase  is  most  nearly  stationary  may  he  chosen  as  origin.  In  any  event,  the  ex- 
pression for  the  phase  in  equation  (62)  can  he  considerably  simplified  hy  a 
translation  and  rotation  of  axes  in  the  xy-plane.  Thus  we  introduce  the  new 
variables  x'  andy',  defined  by  the  equations: 

x'  =  X  cosy  +  y  sin^  +  H,   y«  =  -x  sin)<'  +  y  cosj^  +  K, 

where  tan  2yi    =   = —  ,  and  H  and  K  are  so  chosen  as  to  eliminate  the  first- 

degree  terms  in  the  transformed  equation. 

If  we  set 

A'  =  A  cos^^   +  B  cosV'  sin)^   +  C  sin^^ 

C  =  A  sin^^   -  B  cosji'   sin  y  +  C  cos^J^ 

Ml  =  A'H  +  C'k^  -  D'H  -  B'k  +  M 

where  D'  =  D  cosy   +  E  siny*   ,  !•  =  E  cos)^  •-  D  sin^'  , 

D'  E'  /  /■    \ 

then  H  =  -rr-f  ,   K  =  -^jrr  •  With  these  substitutions,  equation  (62;  becomes 

I(.^.r^)  =/f(r.).^'^'-''*='^''*«>d..dy.         (63) 

It  can  be  seen  that  the  origin  of  the  x',  y'  coordinate  system  is  a 
point  of  stationary"  phase  for  the  integrand  of  equation  (63).   If  there  exists  a 
stationary  phase  point  on  the  surface,  in  accordance  with  earlier  remarks,  it  would 
have  been  chosen  as  the  origin  of  x,y  coordinates  and  hence  also  of  the  x',y' 

coordins.te:=,  the  change  of  variables  consisting  merely  of  a  rotation  in  this  case. 
On  the  other  hand,  if  ther-  is  no  stationary  phase  point  on  the  surface,  the  origin 
of  the  x,y  coordinates  is  a  non-stationary  point,  and  the  origin  of  the  x',y' 
coordinates  is  then  in  the  plane  tangent  to  the  surface  at  this  point.  In  this 
case  the  origin  of  the  x',  y«  coordinates  will  normally  lie  outside  the  range  of 
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integration  in  equation  (63),  and  there  will  then  be  no  stationary  point  in  the 
region  of  integration.   It  is  therefore  to  he  expected  that  the  approximate  value 
to  he  obtained  for  the  integral  —  which  is  small  in  this  case  —  will  he  less 
accurate  than  that  ohtained  when  there  is  a  stationary  point  on  the  surface. 

When  there  is  a  stationsiry  point  on  S  it  is  to  he  expected  that  the 
integral  will  he  the  greater  the  more  nearly  constant  is  the  phase  in  its  neigh- 
horhood.  Thus  if  either  A'  or  C,  the  second  derivatives  of  the  phase  in  equation 
(63),  vanish,  one  would  expect  that  the  integral  will  he  large  in  absolute  value, 
and  still  larger  when  both  vanish.   The  vanishing  of  these  quantities  depends  upon 
the  positions  of  r.  and  r  and  on  the  curvatures  of  the  surface  at  the  stationary 

point,  as  can  be  seen  from  their  definitions.  When  these  quantities  are  such  that 
neither  A'  or  C'  is  zero,  the  following  trajisformation  can  be  made: 


Equation  (63)  then  becomes 


=  yoT,, 


K 


dy"     . 


(6U) 


(65) 


where  J  is  the  Jacobian  of  the  transformation,    and  is  given  by 
.  -     3(x".y")     _  2   L...|4     -     2     L^         B^|2 


2 
n 


1,1  1  ^  .   .  22,1 

D^^2  ~  ""^l  "  2  ^^       5^2 


o  -  az. 


^ 


1^2'UVi       ^         l^lV^g, 
-2r2^J[^2-^^-2^>B^ 


2'      bz-'^ 


Let  kjK^=K^,>  5-    =   sinJr   cos(0       ,   — 


1 
sin  K"   sint^      ,   5-    =  cos  JT      j 


2^ 


(65') 


2  '?  2 

^    =  sin  «*.     cos   0,       z=r-    =  sin  ok  sin  0   ,     =^  =  cos «.    ;    a  +  b 

2  2  ^2 

:c.2  y  2 

l^ab  =  G     ;    2a  -i^    +  2b  ^    =  G     (9)   sin^iT   . 
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Then  J  has  the  form  given  in  equation  (U^)-),  page  24. 
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Prom  equation   (65)  aucL  the  first  expression  for  J  it   is    q)parent   tJat 
the  integral  is  large  v;hen  either  or  "both  of  A'   and  C   are  small.     For  values  of 
r     such  that  A'  or  C    or  toth  sure  zero,   the  integral   in  eqaatlon  (65)  "becomes  infin- 
ite.    An  investigation  of  such  points,   given  in  section  5i    shows  that  they  are  the 
positions  of  conjugate  points  when  only  one  of  A'   and  C*  vanishes,   and  point  taa.gee 
or  foci  when  hoth  vanish.     Actually,    equation  (65)  is  xiot  valid  in  such  cases,   "but 
it  nevertheless  indicates  that  the  integral  hecomes  large  in  absolute  value  near 
these  points.     Equation  (63),  which  must  he  used  at    such  points,   shows  that  the 
integral  remains  finite  there.     The  simultaneoas  vanishing  of  A'    and  C,  which  is 
necessary  for  the  existence  of  point   images  or  foci,   occurs  when  certain  relations 
involving  the  curvatures   of  the  surface,   the  distances  of  the  source  and  ohservei- 
tion  points,   and  the  angle  of  incidence,   are  satisfied.     At  these  points  the  'ields 
are  large.     Point  images   or   foci  at  which  the  fields  are  still  greater  may  "be  ohtained 
if  higher  derivatives  of  the  phase,   say  third  derivatives,  vanish.      Vanishing  of 
these  derivatives  will  require  the  surfewie  and  the  incident  ray  to   satisfy  additional 
conditions  v/hich  have  not  yet  heen  investigated. 

In  order  to  evaluate  the  integral  in  equation  (65)  the  limits  of   Integra^ 
tion  mast  he  determined.      The  region  over  vdiich  the  integral  extends  is  the  pro- 
jection of  the  surface  S   on  the  x",   y "-plane.      If  this  region  is  a  rectaiigle  with 
sides  parallel  to   the  x",  y"-axes,   or  a  circle  with  center  at   the  origin,   the 
integral  can  he  evahiated  explicitly  provided  the   slowly  varying  quantities  f(r') 
and  J  are  replaced  hy  their  average  values.     Any  other  region  will  he  approximated 
hy  a  rectajigular  region.     This  approximation  might  he  expected  to  yield  good  results, 
since  it  agrees  with  a  more  accurate   eveilxxation  if  the  region  is  small,   when  the 
integral  depends  only  upon  its  area  and  not  on  its  shape;     and  if   it    is  large,   when 
the  integral  is  usually  independent   of   the  shape.     Thus,   the  limits  of   integration 
will  he,    in  the  rectsuagular  approximation,* 


V  =  — -    <  ^  <  i;  = 

0      I — — -     -      -    1 


(66) 


Vo  = ^     7<\=  ^' 


rir  "■       rir 


* 

If  the  houndary  is  a  circle  of  radius/*,  in  the  x'',y"  coordinates,  with  center  at  the 
origin,  and  A'/\A'l  =  C'/|C'|   =  +  1,  the  integral  yields 

lir^,r,)^L^     .  2i(l-e-^^>''/2) 


Here  »q,   ^1,^01  ^1   ^^®  ^^   corresponding  limits  in  the  x',  y'  coordinates, 
and  have  the  dimensions  of  length. 
Equation  (65)  now  becomes 

/-U  +   TT  ,  2 

where  r.(u)  is  the  Fresnel  integral  F   (u)  =./     e       ?  d7. 

In  the  firwt  hrsicket  of  equation  (67)  the  sign  is  chosen  the  same   as  that  of .  7-ji  , 

C 

and  in  the  second  tracket  the  same  as  that  of  -rprn-  .  This  is  the  approximate 

value  of  the  integral  in  equation  (U2)  which  will  he  applied  to  the  approximate 
©valuation  of  the  integral  for  E/, \(r)  in  equation  (Ul). 

From  the  expression  for  iCr^^.r  )  in  equation  (67)  the  following  conclusioas 

can  he  drawn.  If  the  rectangle  is  symmetric  ahout  the  origin  in  the  x,y  coordinates, 
with  length  L  and  width  W, 

then 

In  this  case,  if  variations  in  the  slowly  changing  quantity  J  are  neglected, 

1{t^,t   )       has  a  maximum  with  respect  to  r^  or  r  when  H  =  K  =  0,  which  occurs 
if,  and  only  if,  D  =  E  =  0.  From  the  definitions  of  D  and  S  this  means 

iVi  ^  V2  Vi  _    V2 

Equations   (69)  are  equivalent  to  the  law  of  reflection  if  It,   =  k   ,    and  to     Snell's 
law  if  kj^  ^  k^.     In  directions   other  than  those  given  ttg-  these  laws,   either   or 
hoth  of  H  and  K  are  different  from  zero.     Since  they   are  linear  expressions   in  fc, 
and  k^,   they  hecome  infinite  vri. th  k^  and  k  .     Then  from  equations    (68),   either 

»o  ^^   *  1»   °^    I0  ^^  "^1'   °^  "^o^^,   hecome  infinite,   and  the  difference  of  at 
leasb  one   of   the  pairs   of  Fresnel  integrals  becomes   zero.     Thus,    in  the   limit   of 
zero  wave  length,    in  which  case  k^  =  k^  =  00    ,    the  integral   is  zero  except  in  the 

directions  given  hy  geometrical  optics. 
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Appendix  II 


In  this  appendix  the  results   of  the  present  paper  for  the  reflected  and 
transmitted  fields  v/ill  te  specialized  to  the  case  of  a  plane  wave  incident  on  an 
infinite  plane  plate.      They  will  then  he  compEired  with  the  eawict   solution  of  this 
prohlem  given  hy  R.   K.    Luneherg. 

Since,    in  the   case  of  an  infinite  plane,    every  point  receives  either 
specular  reflection  or  direct  transmission,   one  has  ««.=»'    ,  Q  =f    for  reflection 
and  «*  =  TT-r     ,  0  =  IP    for  transmission.     The  curvatures  G     and  G     are  "both  zero 

at  every  point  on  the  surface  and  therefore  the  directions  u_    ajid  u     are  arhitrary. 

For  convenience,  u,    will  he  taken  in  the  plane  of  incidence,  positive  in  the 

direction  of  the  reflected  ray.  By  definition,    the  n-direction  is  normal  to   that 

surface   of  the  plate  nearest   the  source,   and  is  positive  into  the  plate.     The 

■n„-direction  is   then  determined, 
2 

For  simplicity,   an  incident  field  with  the   electric  vector  normal  to  the 
plane  of  incidence  vd''-l  he  assumed.      The  case  in  which  the  electric  vector  is 
parallel   to  the  plane   of  incidence  may  be  treated  in  exactly  the   same  way  "by 
noticing  that  in  that   case  the  magnetic  vector  is  normal  to   the  plane  of  incidence. 
The  following  calculations  and  results  then  apply  provided  one  interchanges  B  with 
H  and  ('  with  -/* .     Thus   a  plane  wave  of  arhitrary  polarization  may  he  dealt  with 
since  it  can  always  he  )^r^itten  as  the  sum  of  two  fields,    one  with  a  normal   electric 
vector  and  the  other  vith  a  parallel  electric  vector. 

A  plane  wave  with  its  electric  vector  normal  to  the  plane  of  incidence  may 
he  written  ^^   ,^_^  , 


a    o 

2    0 


=  A.  e         ,  E^  =  B^  =  0 


#-7T—       ik„  I  r-r  t 
ti-»/a   .y-    a    o' 

H  -  A  J smo  e         , 

n    0  '  14  * 

/^  a 

rr-  ik  Ir-r  \ 

H,  =  -A  7—^  cos  y  e  ^    °   ,  H„  =  0 
1      0'  ly  '2 

Here  A  is  a  constant  scalar,  the  amplitude  of  the  field.  The  source  is  located  at 
r  ,  which  recedes  to  infinity  along  a  line  in  the  n,  u, -plane  making  the  angle  i^ 

with  the  n-axis.   Let  the  coordinates  of  r  he  z  ,  x  ,  0  in  the  n.xL,  ,u_-directions. 

0  00  '     A.       2 

respectively.     Those   of  any  other  point  r  may  he  written  as  z,x,0,   since  the   field 
and  reflecting  ohject   do  not  vpjy  v/ith  Up  and   therefore  any  point  may  he   taken   to 


^. 


"be  in  the  plane  u,  =  0.  One  then  has  for  the  distance  |r  -  r  | 


1/2 


|r  -  rj=   U^-^o^  -*■  U  -   X  )^]    =  (z  -  z  )  cosr  +  (x  -  x  )  alnr 


The  phase  ^  of  the  plane  wave  is  so -adjusted  that,  as  r  goes  to  infinity, 

?  -  k  z  cosiT  -  k  X  slnf  is  an  even  multiple  of  TT.  One  then  ohtains  the 
^  °         ^     °  ik  (z  cosr  +  X  sini  ) 

usual  expression  for  a  plane  wave;   e  .  We  oculd,  of  cexirse, 

have  assumed  this  at  the  outset;   the  foregoing  discussion  justifies  this  form, 

Trom  the  geometry,  one  finds  that  ^  ~  p  -  <"   and  «  =  ^  . 

Using  the  above  data  in  equation  (^5)  ^^^   *^  corresponding  equation  for 

H.  one  finds  for  both  the  reflected  and  transmitted  fields  at  any  point .that 
•  Ik  (z  cosr  +  X  slnir  ) 

E,  =  B  =  H„  =  0.  The  phase  factor  becomes   i  e  ^  for 

ik  (-Z  cost"  +  X  sin  I'  ) 
transmission  and  1  e  for  teflection,  since  in  both  cases 

A  and  C  are  positive  and  the  Fresnel  integrals  have  their  asymptotic  values. 


The  geometrical  factor  becomes 


^1*^2 


=  -  1 


for  transmission,  and 


^* 


=  1  for  reflection,  since  D.,  Is  infinite. 


Using  these  expressions  in  equqtion  (U5)  and  the  corresponding  equation 
for  H,  the  remaining  reflected  and  transmitted  field  components  are 
Reflected         _^^ 
So 
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The  ex£u:t   solution  of  reflection  and  transmission  of  plane  waves  "by  an 
infinite  plate  of  any  thickness  has  "been  given  "by  S.   K.  Luneherg       (ll)  p.^^      , 
The  amplitude  of  the  reflected  field  is   that   of  the  incident  field  multiplied  "by 
the  reflection  coefficient  R;    similarly  for  the  amplitude  of  the   transmitted  field 
and  tlE  traJismission  coefficient  D.      TEnese  coefficients  are  given  helow  for  the 

case  of  a  parallel  electric  field.      The  sign  of  i  has  "been  changed  in  order  to 

-i«^t 

yeild  agreement  with  our  result  since  we  use  the  factor  e      and  Luneberg 

+114'  t 

has  used  e  .     Then  one  has 
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The  approximation  is  valid  when  M.T     is  small,  which  is  the  case  if 

h  is  small  enougla.     By  reference  to  articles  k  and  5  of  Luneherg's  p^er     one 

sees  that  the  dependence  of  the  reflected  and  transmitted  fields   on  x  y  z  is 

the  same   as  that   obtained  in  our  solution.     However,   Laneoerg's   transmitted 

ikg^  h  cos  Y 
field  has   an  extra  factor  e  .     In  order  to   compsire  the  solutions,    this 

must  he  expanded,  yielding 


D  •  =  D  -ik  h  cos  V'  =  1  + 
P  pa 


ik 


"    1  K<^' 

a     D 


<^aH 


cos 


r+ 


^hA-</^a 


By  comparing  these  solutions  with  the  results  of  this  paper,    it  is  seen 
that   the  agreement   is   complete  to  the  order  of  approximation  retained.    It  must 
"be  noticed  that  Luneherg's  result  for  the  parallel  electric  field  has  "been  given; 
to  compare  it  with  our  result  for  the  normal   electric  field  f •    and  -yU  mast  he 

interchanged. 
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List  of  Symbols 

Gaussian  units  are  used  throughout 

2 
^   =  Laplacian  operator 

k  =  — £/t 

c 


14?  =  angular  frequency  of  field 

^  =  dielectric  constant 
/^   =   magnetic  permeatility 

E  =  electric  field  (vector) 

H  =  magnetic  field  (vector) 

c  =  velocity  of  light  in  vacuum 

r,r'  =  position  vectors 

G(r,r')  =  Green's  function  defined  on  page  6 

|r  -  r'l   =  distance  from  r  to  r' 

<J  =  any  function  of  r,  defined  on  page  6 

V  =  three  dimensional  region,  defined  on  p.  6 
S  =  surface  bounding  V 

n  =  normal  vector  to  S  into  Y 

E  ,H  =  electric  and  magnetic  fields  in  absence  of  obstacle,  i.e.  incident  fields 

0   0 

V  =  three  dimensional  region  in  region  (a) 
S  =  surface  of  V 

S'  =  obstacle  surface;  i.e.  surface  of  region  (b) 

(J-  =  sphere  of  very  large  radius  with  center  in  region  (b) 
(J~    =   small  sphere  vfith  center  at  r 

E  ,E  =  limits  of  B(r)  as  r  approaches  the  surface  S*  from  media  (a)  and  (b) 
respectively 

H  ,H  =  limits  of  fi(r)  as  r  approaches  the  surface  S'  from  media  (a)  and  (b)  respectivell 
n,  as  subscript,  indicates  normal  component  of  a  vector  at  the  surface  S' 
t  indicates  tangential  component  of  a  vector  at  the  surface  S' 

u.  =  parameter  which  is  constant  along  one  family  of  lines  of  curvature  on  S' 
Ujj  =  parameter  which  is  constant  along  one  family  of  lines  of  curvature  on  S' 
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x(  u.  Up)  =  vector  function  of  u,  ,  u,  giving  the  surface  S* 

g^ ,  g^  ~     metrical  coefficients  in  element  of  axe  length  (page  10) 

i, ,  ip  =  unit  vectors  tangent  to  lines  of  curvature  on  S' 

E,  ,E-  =  components  of  E  in  directions  of  iitip  respectively 

E.  ,H  =  components  of  H  in  directions  of  i-,,ip  respectively 

G-  =  mean  curvature  at  a  point  of  S'  (page  12) 

S-,  Sjj  =  paradlel  surfaces  of  the  olDstacle  ;  parts  of  S' 

r_,  r'  =  points  on  S^,  S-j.  respectively 

A(r»)  =  ab"breviation  for  integrand  of  surface  integral  in  eq.  (27)  or  (28) 
h  =  thickness  of  obstacle 

y)  -   variable  of  integration  ranging  from  0  to  h 
n'  =  ujiit  normal  to  St  at  r«,  into  obstacle 

•'  -  summation  index,  ranging  from  0  to  00  thru  positive  integers 
E/^\(r),  H/  v(r)  =  coefficients  of  h   in  series  expansion  of  E(r)  and  H.(r) 

in  medltun  (a).  They  are  vector  functions  of  r. 
C(/f)^^tOf)=   vector  aatplitudes  of  incident  electric  and  magnetic  fields 

^  t  P  t  ^   ,  ^     -   angles  defined  on  p.  21 

f  '^  j" ^  V  ^"  ~  limits  of  integration  defined  on  p.  23  and  p.  U3 
F^(x)  =  Fresnel  integral  defined  on  p.  23 

G  =  Seometrlcal  factor 
P  =  Phase  faxjtor 

The  symbols  in  the  following  group  are  defined  on  page  2k. 

M« 

^1 

^2 

J 

0,4) 

G 
g 

G 

Several  symbols  not  spearing  in  the  above  list  are  defined  in  the 
appendices  and  used  mainly  in  them. 
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